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ABSTRACT 

We have conducted a detailed object-by-object study of a mass-complete (M* ^ 
10 11 Mq) sample of 56 galaxies at 1.4 ^ ^ 2 in the GOODS-South field, show¬ 

ing that an accurate de-blending in MIPS/24um images is essential to properly assign 
to each galaxy its own star formation rate (SFR), whereas an automatic procedure 
often fails. This applies especially to galaxies with SFRs below the Main Sequence 
(MS) value, which may be in their quenching phase. After that, the sample splits 
evenly between galaxies forming stars within a factor of 4 of the MS rate (~ 45%), 
and sub-MS galaxies with SFRs ~ 10 — 1000 times smaller (~ 55%). We did not find 
a well defined class of intermediate, transient objects below the MS, suggesting that 
the conversion of a massive MS galaxy into a quenched remnant may take a relatively 
short time (< 1 Gyr), though a larger sample should be analyzed in the same way 
to set precise limits on the quenching timescale. X-ray detected AGNs represent a 
~ 30% fraction of the sample, and are found among both star-forming and quenched 
galaxies. The morphological analysis revealed that ~ 50% of our massive objects are 
bulge-dominated, and almost all MS galaxies with a relevant bulge component host an 
AGN. We also found sub-MS SFRs in many bulge-dominated systems, providing sup¬ 
port to the notion that bulge growth, AGN activity and quenching of star formation 
are closely related to each other. 

Key words: galaxies: evolution - galaxies: high-redshift - galaxies: structure - galax¬ 
ies: formation 


1 INTRODUCTION 

The transition from a star forming regime to passive evo¬ 
lution, that we call quenching, is perhaps the most impor¬ 
tant event that may happen to a galaxy in the course of 
its life. Yet, the physical processes causing quenching are 
still poorly understood, making this a main frontier issue 


sive) 

galaxies appear at redshifts ~ 2 — 3 (e.g., 

Cimatti et al. 

2004 

Williams et al. 20091: Gobat et al.ll2012l: 

Cassata et al. 

2013 

and that their fraction and number density increase 


steadily with cosmic time, to the point that they come to 
domi nate the ste lla r ma ss in galaxies in the local Universe 
(e.g., iBaldrv et al.l [20041 ). We also know that the quenched 
fraction increases steeply with galaxy mass independently of 
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local environment, as well as with local environmental over¬ 
density, independently of their mass. We therefore speak of 
mass quenching and environment quenching as two distinct, 
separable processes dPeng et al ] 201 OhJ . 1 201 2jf . Yet, the phys¬ 
ical nature of such processes is still merely conjectural and 
it has also been suspected that they coul d be different man¬ 
ifestation of a same underlying process dCarollo et al. I [20141 : 
iKnobel et al.ll2015l ). 

The probability (or relative frequency) of mass quench¬ 
ing appears to increase exponentially with stellar mass 
(M*), or almost indistinguishably with the star formation 
rate (SFR), given the tight SFR-M* correlation that ex- 


ists for main sequence (MS) galaxies (Daddi et al. 

to 

o 

o 

KJ 

Elbaz et al. 20071: 

Noeske et al.l 2007i: Schreiber et al. 

20 ll) 


A variety of radically different options are currently enter¬ 
tained for the mass quenching process, whereby quench¬ 
ing is either an internal or an external process. In one 
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option for the former case sudden energy/momentum re¬ 
lease from star formation and/or AGN (feedback) results 
in the ejection of all gas from galaxies that then turn pas¬ 
sive^ tlTe_) 2 uasarmode’ quenching in current jargon (e.g., 
iGranato et al.ll2004f l. Powerful AGN jets may also heat the 
circumgalactic medium to high temperature thus prevent¬ 
ing further acc retion of cold gas, t he so-called ‘radio mode’ 
AGN feedback (ICroton et al.ll2006h . In another option for an 
external process, the circumgalactic gas is shock-heated to 
high temperatures as the mass of the host dark matter halo 
(Mh) exceeds a critical threshold (of order of ~ 10 12 Mq), 
and therefore it stops to cool and flow in to the galaxy, thus 
disco ntinuing to feed star formation le.g.. lBirnboim fe Dekell 
l2003h . Finally, the growth of a central mass concentration 
(bulge) may quench itself, with increasing shear (differential 
rotation) suppressing the disk instability to form actively 
star-forming clump s, the so-called gravitational (or morpho - 
logical) quenching dMartig et al.|[200^ : iGenzel et al.ll2014al l. 
So, we have at least four options for the physical nature 
of mass quenching. Moreover, the tight correlations existing 
between halo mass, stellar mass, bulge mass and mass of 
the central black hole make difficult to disentangle between 
these options, as all resul t in a correlation of the quenched 
fraction with galaxy mass (IPeng et al.ll2012h . In other words, 
it is not clear whether the pertinent mass in mass q uenching 
is M * or Mh, or the mass of the central bulge dLang et al.l 
l2014h , or even the mass of the central black hole. 


The frequency of AGN activity is found to in- 
crease strongly wi t h galaxy mass, especially at z > 1 


dReddv et al 

TT2005S: Daddi et alj 2007a 

: iFiore et al.l l2008i: 

Brusa et alj 

20091: Boneiorno et al. 

2012: Mullanev et al. 

2012a.b: Cimatti et al. 2013: Forster Schreiber et al. 20141: 

Genzel et al. 

2014bl: Perna et al.1 20151: Rodiehiero et al.l 


20151 '). which provides circumstantial support for a connec¬ 


tion between AGN feedback and mass quenching. However, 
finding the smoking gun of AGN quenching proved to be 
very difficult, as most AGNs are hosted by actively star 
forming galaxies showing no signs of being quenched. Actu¬ 
ally, the mere galaxy-AGN co-evolution concept argues for 
both growing and be quenched together. On the other hand, 
supermassive black holes do not grow in one shot. Rather, 
a galaxy may experience ma ny nuclear activity cycles inter¬ 
leaved with inactiv e phases llNovak, Ostriker fc Ciottal2012l : 
iHickox et alj|2014h . with the former ones leading to feed¬ 
back and driving nuclear outflows, but failing to quench. 
Hence, if AGN quenching really works, not every AGN cy¬ 
cle leads to quenching but just a final, fatal one. Later, once 
star formation is quenched by ejecting the gas, so is the 
nuclear activity as well. This AGN variability on a variety 
of timescales clearly makes the conn ection AGN-quench ing 
difficult to unambiguously recognize dHickox et al.|[20l3) . 

While mass quenching works equally for satellites 
and centrals, environment quenching appea rs to be exclu- 
sive of satell i te gal axies dPeng et al.l 120121 '). but see also 
iKnobel et al.l |2015). Also for the physics of environment 
quenching vario us options exist. Ram -pressure stripping is 
one possibility dGunn fc Gotti 1 19721 '). perhaps favored by 
the finding that the quenched fraction of satellites corre¬ 
lates better with the local overdensity inside groups than 
with group richness, a proxy for the mass of the host 
halo dPeng et al.l 120129 . But other options include stran¬ 


gulation dLarso n. Tinsley fe Caldwelllll980l ') and harassment 
dMoore et al. 19961 '). 

In the Peng et al. phenomenological model, mass 
quenching starts promptly at relatively high redshifts and, of 
course, acts preferentially on massive galaxies. Environment 
quenching sets in progressively, becomes more important at 
lower redshifts, following the growth of overdensities and 
large scale structures. 

In this paper we focus on mass quenching, and to do so 
we deal only with the most massive galaxies, namely those 
with M* ^ 10 11 Mq. We then explore the relatively nar¬ 
row redshift range 1.4 ^ z ^ 2, because this corresponds to 
an epoch when mass quenching has started to work at full 
steam, while environmental quenching is just about to be- 
gir0- All such massive star forming galaxies must soon start 
to be quenched, otherwise -keeping to form stars at the main 
sequence rate - would soon r esult in a dramati c overgrowth 
of their mass dRenzinill2009l : IPeng et al1l2010bh . So, all such 
galaxies are almost immediate precursors to quenched, pas¬ 
sively evolving, red and early-type galaxies. Thus, a mass 
selection M* > 10 11 Mq of galaxies at the mentioned red¬ 
shifts ensures that it will include only galaxies which are 
either already quenched, or being in the course of quench¬ 
ing, or that will soon (i.e., < 10 9 years) enter the quenching 
phase. 

Such mass and redshift range has been already 
widely explored by many observational studies using large 
databases, in par ticular tryin g to a ddress the quench¬ 
ing issue (e.g., Wuvts et al. ] l 201 ll: iBarro et al.l I2OI3I : 

I Williams et al.ll2014f) . Some of these studies deal with very 
large numbers of galaxies to ensure the statistically signifi¬ 
cant samples that are crucial to understand the global trends 
of galaxy evolution. Given the sheer size of the samples, 
automated methods are used to derive a few fundamental 
parameters for each galaxy, such as SFR, mass, structural 
parameters etc. As such, not all information present in the 
data could be used. We then restrict ourselves to work, in¬ 
tensively rather than extensively, on a quite small sample 
of galaxies, examining them one by one. The purpose is 
first to unambiguously distinguish between galaxies which 
are forming stars at the main sequence rate, those which 
are already quenched, and those which appear to be still 
star-forming, but at a rate significantly below the main se¬ 
quence, which may be caught in the quenching phase. We 
also aim to investigate possible links between galaxy specific 
star formation rate (sSFR), morphology, AGN activity. The 
presence of a we ll known overdensity a t 2 ~ 1.61 i n the 
GOODS-S field (ICastellano et ahl 120071 ; iKurk et al.l 120091 : 
IGiavalisco et al.M201lh also allows us to check the relative 
abundance of quenched galaxies, bulge-dominated galaxies, 
and AGNs within and outside the overdensity. 

The paper is organized as follows. Section [2] describes 
the sample selection and the used catalogs and data-sets. 
In Sections [3] we derive galaxy stellar mass and SFR and 
define the M* —SFR main sequence at 1.4 ^ z ^ 2, used 
as reference throughout all the paper. AGN hosts are dis- 


1 For example, in the case of the lPeng et al.l j~2Q I Oil ) phenomeno¬ 
logical model at z ~ 2 and M* = 10 11 Mq the mass quenching 
rate is over two orders of magnitude higher than the environment 
quenching rate, see their Figure 13. 
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cussed in Section El In Section E galaxies are classified in 
quenched or star-forming, based on their position relative 
to the main sequence, multi-wavelength information, col¬ 
ors, and SED-fitting analysis. In Section Ewe perform mor¬ 
phological analysis based on single- and double-component 
surface brightness fitting, and separately treat some pecu¬ 
liar system. Based on our results, we discuss the role of 
the AGN activity, and of the environment, in quenching 
star-formation in Sections E and [7] respectively. Section E 
presents the conclusions. Throughout the paper we assume a 
ACDM cosmology with Ho = 70 km s -1 Mpc^ 1 , Am = 0.27, 
and Ha = 0.73. All stellar masses and SFR are quoted for a 
ISalpeteil dl955l l initial mass function (IMF), and magnitudes 
are given in the AB photometric system, unless explicitly 
stated otherwise. 


2 DATA AND SAMPLE SELECTION 


The relatively small field (10 / x 16') of the Great 
Observatories Origin s Deep Survey-South (GOODS-S, 
iGiavalisco et al.l 120041 1 represents the optimal choice for 
studying in detail a representative sample of the most 
massive galaxies at 1.4 Sj z ^ 2. In fact, it pro¬ 

vides at the same time both high-resoluti on images in 
the near-IR with WFC3 /HS T (CANDELS, IGrogin et all 
l201ll ; iKoekemoer et ahl 201l|) , an d the deepest data in 
the mid-IR (S vitzer. iRieke et al.l 12 00411. far-IR ( Herschel. 
PEP/HerMES, lElbaz et al .1120111: iLutz et al.ll201lll. and X- 
ray b and (Chandra Deep Field-South, CDF-S. Ixue et ahl 
l201ll l. The corresponding database is crucial to study the 
optical rest-frame morphology, accurately estimate the SFR 
for the most massive objects (most of which are highly 
obscured by dust), and to unveil X-ray excess due to the 
AGN contribution. Last but not least, the large amount 
of ground-based (VLT) and space-based (HST, Chandra, 
Spitzer) data from the GOODS Treasury Program, and the 
availability of spec troscopic redshifts from the GOODS and 
GMAS S progr ams (Mm^ejka^talJgOO^ ICimatti et ai1l2008l : 


IPopesso et al.ll200£ : Kurk et al.l 20131 ) enable us to charac¬ 

terize the studied galaxies in an optimal way as currently 
possible. Of course, the resulting sample will be too small 
for coping with statistics or cosmic variance, but this is not 
a limitation in our intensive approach. 

I 11 this work we focus on the 56 most massive (M* ^ 
10 n MQ) galaxies at 1.4 ^ z ^ 2, 31 of which have spectro¬ 
scopic redshifts. Two more objects entered the original sam¬ 
ple, but they have been excluded from this study, because 
their WFC 3/HST images are not available, or too noisy to 
perform surface brightness fitting, due to the proximity of 
saturated stars. The sample was culled from the A- s elected 
(A(Vega) < 22) multi-band catalog of lDaddi et al.l (l2007bl . 
hereafter, D07), including data from all the available filters 
in GOODS-S, i.e., HST /ACS optical, F435W (B), F606W 
(V), F775W (I), and F850LP(z), VLT/ISAAC near-IR, J, 
H, K, and Spitzer/IRAC, 3. 6, 4.5, 5.8, and 8.0 qm (for more 
details on the data-sets see IGiavalisco et al.ll2004l 'l. The op¬ 
tical/near IR photometry was then complemented with the 
24 p m catalog (Daddi et al., in prep.), built as summarized 
in Se ction ED from the Spitzer/MIPS images dR.ieke et al.l 
120041 1. Far-IR fluxes were extracted (Daddi et al., in prep.) 
from the publicly released PACS 70-160 pm data from Her¬ 


schel GOODS (H-GOODS, lElbaz et al.ll201lll , and SPIRE 
250 pm data f rom Herschel Multi-tiered Extragalactic Sur¬ 
vey (HerMES, lOliver et al M3)- For objects without spec¬ 
troscopic information we used photometric redshifts from 
the public GOODS-MUSIC catalog, which agree well with 
the spectr oscopic ones (Az/(1 + z) ~ 0.03) for galaxies at 
z < 2 (cf. iGrazian et al.l I20061 . 120071 1. Although the uncer¬ 
tainties on photometric redshifts could result in the inclu¬ 
sion in the sample of a few lower-redshift contaminants, we 
avoided applying color criteria to pre-select high-z objects, 
in favor of completeness. The AST/WFC3/F160W 17-band 
image mosaic (> 5<r point source sensitivity for H 160 < 27.7, 
AB system) drizzled to a pixel scale of 0V06, was exploited 
to study galaxy morphology in the optical-rest frame with 
a very high resolution (FWHM~ 0'/l8 ~ 1.5 kpc). For 
more details on th e ob servations and da t a red uction see 
IGrogin et al.l (1201 il l and IKoekemoer et al.l d201 ill . 


2.1 MIPS fluxes 

The MIPS source extraction on the whole GOODS-S 
MIPS/24 pm image mosaic will be described in a future 
paper (Daddi et al., in prep.). Here we briefly summarize 
the main issues, useful for the purpose of this work. The 
24 /rm/MIPS counterparts were identified, based on the 
IRAC (3.6 pm) prior posi tions, usi ng a PSF fitting method, 
with GALFIT (version 3, IPeng et al.ll2010al l. The measure¬ 
ment errors were calibrated on direct Monte Carlo simula¬ 
tions, inserting one artificial source at a time in the real 
image, and measuring it together with all other priors. The 
catalog was built in such a way that a MIPS flux, or a flux 
upper limit, in case of MIPS-undetection, is associated to 
each IRAC source. The use of IRAC priors helps with the 
deconfusion of objects blended in MIPS, but is still biased 
against sources which are also blended in IRAC. In the D07 
catalog, IRAC sources are associated to the A-band coor¬ 
dinates using an encircling radius of 270. However, when 
the angular separation between the IRAC and A-band po¬ 
sitions is A6k-irac > 075, the galaxy is probably blended 
in IRAC (cf., D07), and generally disregarded. Since ob¬ 
jects with A Qk-irac > 0”5 represent a not negligible frac¬ 
tion in our sample (i.e., 6/56), we used prior positions from 
AS , T/WFC3/F160W to deblend them, as described in Sec¬ 
tion [4] and Appendix 1X1 

We considered as MIPS-detection all sources with S/N^: 
3 (i.e., detected above 3<r). For undetected sources, we de¬ 
rived upper limits from the flux errors at 2a or 3cr, for 
sources with S/N< 1 and 1 ^S/N< 3, respectively. In Fig¬ 
ure [p we compare our r esults with those from the catalog 
of iMagnelli et all (120091 . hereafter M09). The figure shows 
the ratio between the MIPS fluxes from our catalog and 
those from the M09 catalog (F24/F24(M09)) as a function of 
F24(M09). For sources undetected in M09 we estimated 2a 
or 3 ct flux upper limits as we did for objects in our catalog. 
With one exception, all the sources detected in our catalog 
are also detected in M09 and are shown as red filled circles in 
Figure E For sources undetected in both catalogs the black 
arrows show the corresponding upper limits in F24(M09) 
and the vertical axis show the ratio of the upper limits in 
the two catalogs. The exception, object #5530, is shown as 
a red arrow. The red (black) shaded histogram show the 
corresponding distributions of flux ratios (of the ratio of up- 
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Figure 1. The ratio between our 24 /rm/MIPS fluxes (F24, Daddi 
et al. in prep.) and those from M09 (F24(M09)), as a function 
of F24(M09). For undetected sources (S/N< 3) we plot 2 or 30- 
flux upper limits, as detailed in the text. Objects detected (un¬ 
detected) in both catalogs are shown as red filled circles (black 
arrows). The object #5530 detected in our catalog, but unde¬ 
tected in M09, is shown as a red arrow. The red and black shaded 
histograms (scaled down to arbitrary units), show the distribu¬ 
tions of sources which are detected and undetected in our catalog, 
respectively. Gaussian fits to the histograms are shown as red and 
black dot-dashed lines. 


per limits) for detected (undetected) sources. Gaussian fits 
to the histograms, shown as dot-dashed lines, give a mean 
(fj,) and standard deviation (a) as # ± a = 1.06 ± 0.14, 
and 1.21 ± 0.74, for detected, and undetected sources, re¬ 
spectively. The comparison shows that for sources detected 
above 3a, the 24 /rm/MIPS fluxes are in general agreement 
with those from M09. On average, also the upper limits are 
in fair agreement between the two catalogs, whereas for a 
small number of objects our upper limits are more conser¬ 
vative. 


3 THE MAIN SEQUENCE (MS) OF 

STAR-FORMING GALAXIES AT (z) = 1.7 

It is now well established that a tight relation exists be¬ 
tween the stellar mass (M*) and the star formation rate 
(SFR) for star-forming galaxies , both in th e local Universe 


Elbaz et al. 

200?l. Noeskc et al. 1 . 2007. Pannella et al.l 200fll. 

Karim et al. 

l 201 ll: Rodighiero et al.l 120141: ISchreiber et al.l 


201fll ). At all redshifts the majority of star-forming galax¬ 


ies follow this relation within a < 0.3 dex dispersion with a 
small fraction of outliers (cf. iRodighiero et all 1201 if) , thus 
it is commonly called the main sequence (MS) of star¬ 
forming galaxies. While the slope of the SFR-M* relation 
does not change much with redshift, its zero-point, i.e., the 
specific SFR (sSFR=SFR/M*) at fixed mass gradually in¬ 
creases with redshift at least up to z ~ 2.5 llWuvts et all 
l 201 ll ; I Sargent et al.l [20121 : IWhitaker et al.l| 2012 l ). At higher 
redshifts, it is still debated whether the sSFR remains 
const ant or keeps (sl i ghtly) increasing dGonzalez et al.l 
l 20 ld : IStark et all 120131 : Ide Barros. Schaerer fe Starkl 20141 : 


IBethermin et all 120131 : iGonzalez et al.ll2014l ). So, massive 
galaxies at high redshifts formed stars and grew their 
mass at a much faster rate, fueled by a larger reservoir of 
fresh gas, with resp ect to local galaxies dDaddi et a.1.1 12010 : 
iTacconi et al.ll 2 Qlol) . 

To study our sample galaxies as a function of their posi¬ 
tion on the M*-SFR diagram, we define here the MS of star¬ 
forming galaxies at ( z) ~ 1.7, to be used as a reference in 
quantifying as sSFR/sSFR(MS) the deviation form it of any 
given galaxy. We adopted th e following em pirical relation 
from Bethermin et all d 2012 l ) , based on the ISargent et all 
d2012 ) parametrization of the sSFR evolution with redshift 
for MS galaxies, which in turn is derived by fitting a combi¬ 
nation of mid-IR, UV, far-IR, and radio data: 


sSFRms(z, M * 


sSFRms.o) x 


M* \ Pms 
10 11 M@) 

x(l + Z y MS , ( 1 ) 


where sSFRms.o = 0.063 Gyr -1 is the sSFR at 2 = 0 for 
Af* = 10 11 Mq, /3ms = —0.21 is the s lope of the sSFR- 
M* relation derived bv lRodighiero et al.l d201lh . and 7 ms = 
3 describes the evolution of the normalization of the MS 
out to redshift 2 : = 2. 5. All these paramet ers are g i ven in 
IBethermin et al] d2012l . Table 1) assuming a lSalpeten dl955l ) 
IMF. From Eq. (|TJ) it follows that the MS at (z) « 1.7 is 
given by: 


log SFR = — 6.6 + 0.79 log M*. (2) 


3.1 Stellar Mass 


The stellar masses were computed based on the K-band 
magnitudes, according to the empirical relation derived in 
D04 and D07 (with a typical error of a factor of ~ 2), for 
homogeneity with th e samples used to deriv e the reference 
main sequence (i.e. IRodighiero et aid 1201 ll : ISargent et aid 
120121 : IBethermin et alT 20121 ). in which stellar masses were 
derived using this method. In D04 the empirical relation was 
tested and calibrated on the r esults from best-fit S pectral 
Energy Distributio n s (SE P) of lFontana et al.l (120041 ). using 
iBruzual fe Charlotl (120031 ) stellar population models (here¬ 
after BC03), and Salpeter IMF. As shown in Section 14.21 
we verified that the stellar masses derived with this method 
(M* D04) agree well with those obtained from the SED fitting 
analysis (i.e., M* sed) also for our massive sample. 


3.2 Star Formation Rate 

By combining the deep PACS Herschel data with all the 
high-quality multi-band information at shorter wavelengths, 
we estimated SFRs including both the contribution of the 
UV rest-frame (from the unobscured part of the light), and 
of the total IR emission, i.e., the light reprocessed by dust. 
The total SFR is then indicated as: 

SFR(IR + U V) = SFR(IR) + SFR(U V). (3) 

Following lKennicu ttl (ll998l) . we inferred the UV contribution 
(SFR(UV)), from the 1500 A rest-frame luminosity (T 1500 ), 
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which in turn was derived from the observed B-band mag¬ 
nitude, as done in D07. 

The IR contribution, SFR(IR), was derived from the 
total IR luminosity integrated over the full infrared spec- 
truim i .e., rest - frame L[8-1000 /im] (hereafter LIR) using 
the iKennicuttl (Il998l l conversion, for a Salpeter IMF (i.e., 
SFR(IR) = 1.7 x 10“ lo MQyr _1 Lg 1 ). LIR was measured by 
fitting the Spitzer/ MIPS and Hersc/ieZ/PACS/SPIRE pho¬ 
tometry (i.e., from 24 /im to 250 pm) to t he SEP tem¬ 
plates for MS galaxies at z = 1.325 — 2.25 from lMagdis et all 
d2012f) . The possible contribution of an o bscured AGN was 
also t aken into account in the fit, as in iRodighiero et all 
(1201 ill . The uncertain ties on LIR were estimated from the 
X 2 variations following I Avnil (Il976l l for the case of one single 
parameter. This is a robust method to retrieve the LIR when 
a galaxy is detected in at least two filters (i.e., MIPS/24 pm 
and one Herschel filter). Otherwise, it only provides an up¬ 
per limit to LIR. Hence, the LIR of objects detected in 
MIPS, but undetected in Herschel , was estimated based on 
the 8 pm rest-frame luminosity (L8 = v L^S^un]), accord - 
ing to the empirical relation proposed bv lElbaz et alj (l201ll l 

LIR = L8x4.9t 2 :2-, (4) 

where lcr errors on the coefficient are indicated. For galax¬ 
ies in the considered redshift range (z = 1.4 — 2) the 8 
pm rest-frame luminosity can be derived from the observed 
MIPS/24 pm flux (F24), with a relatively small k-correction 
depending on the redshift (0.8 < k COI r(z) < 1.9), as : 

L8 = V 24 M m F2A^kDL(^z) fccorr W, (5) 

where DL(z) is the Lu minosity distance. T he k-correction 
was derived using the iMagdis et alj (|2012l l templates for 
MS galaxies at 2 = 1.325 — 2.25. For MIPS-detected/PACS- 
undetected objects, we verified that the MIPS-derived LIRs 
are consistent with the 2cr upper limits derived from the IR 
SEDs. 

For objects undetected in MIPS, upper limits to the 
LIR, and correspondingly SFR(IR), were estimated from the 
MIPS/24 pm detection-limit (i.e., accurately measured on 
the MIPS images, as described in Section mu . using Equa¬ 
tions © and ©• We ignore here a possible contribution to 
LIR coming from circumstellar dust in old stellar popula¬ 
tions, which may be relevant for sub-MS galaxies. Hence, 
for such galaxies, the resulting SFR(IR+UV) should be re¬ 
garded as an upper limit, while a better constraint on the 
SFR can be derived from optical/near-IR SED fitting (see 
Section 14.21) . 


3.3 Identification of AGN hosts 

We cross-correlate d our sample w ith both the 4 Ms Chan¬ 
dra X-ray catalog ( |Xue et al.l201lh . and the VLA (1.4 GHz) 
radio catalog dMiller et al. 2008T) . using an encircling radius 
of l'.'O, and found 19 and 5 counterparts, respectively (cf., 
Table [lj . All the radio-detected objects are also X-ray de¬ 
tected, at least in one of the three bands of the 4 Ms Chandra 
catalog (i.e., Soft: 0.5-2 kev, Hard: 2-8 kev, and Full: 0.5- 
8 kev). All the X-ray-detected objects i n our sam p le ha d 
been classified as “AGN candidates” bv IXue et alj (l201lTi . 
However, we considered as AGN hosts only those galaxies 


showing an excess in X-ray and/or radio luminosity with re¬ 
spect to the expected SFR, derived by summing the UV- and 
the IR-derived contrib utions, i.e., SFR(IR+UV). We used 
the lMineo et al.l (120141) calibration, relating the total, inte¬ 
grated X-ray luminosity LX[0.5—8 kev] to the galaxy SFR 
(LX/SFR=(4.0±0.4) x 10 39 [erg s -1 / M 0 yr" 1 ]), and found 
that all the X-ray detected galaxies, but two (i.e., #4705, 
and #5503) show an X-ray excess (at least by a factor of 3). 
In a sim ilar way, by adopting the conversion of llvison et al.l 
toiobiia , we verified that all our radio-detected galaxies 
(#686, #3231, #5503, #5860, and #5898) have a higher 
radio-emission (by a factor > 2.5) than what expected only 
from star-formation, thus revealing the AGN contribution. 
Since #5503 shows radio-excess, but no X-ray excess, we 
consider this object just as an AGN candidate. We notice 
that this source also exhibits a peculiar, asymmetric, mor¬ 
phology, and a large offset (~ 1 arcs) between the optical and 
the X/radio positions, indicating that it could be a merg¬ 
ing system. In the following figures it will be flagged with a 
different symbol (large circle + four-pointed star) with re¬ 
spect to the other more secure AGNs, just marked as four- 
pointed stars. Based on the LX-SFR(IR-fUV) comparison, 
object #4705 is instead considered as a normal star-forming 
galaxy. We conclude that 18 X-ray detected AGNs (17 se¬ 
cure, and one candidate) are included in our massive sam¬ 
ple (see Table [lj. This roughly corresponds to ~ 30% of the 
galaxies with M* ^ 10 11 Mg at 1.4 ^ 2 ^ 2 in GOODS-S. 
Interestingly, we will show in the following that AGN hosts 
are not only among star-forming galaxies, but they represent 
~ 22% of the quenched sources. 


4 SOURCE CLASSIFICATION 

4.1 Position on the M*-SFR plane 

The position occupied by our massive galaxies in the SFR- 
M* plot is shown in Figure© The MS relation at (z) = 1.7, 
derived from Eq is bounded by two dashed lines, repre- 
sentin g offsets in SFR 4x abov e and 4x below the MS fol¬ 
lowing Rodighiero et all (1201 il l, and roughly corresponding 
to the ~ 2.5 a dispersion. Here, and throughout all the paper, 
the sample is color coded based on the galaxy IR properties. 
Blue filled circles are objects detected at Herschel (PACS) 
wavelengths, and hereafter are indicated as the PACS-d sub¬ 
sample. Sources detected at 24 ^im/MIPS, but not detected 
by Herschel , hereafter the MIPS-d galaxies, are shown as 
green filled circles and have, on average, a lower specific 
SFR. (sSFR) compared to the previous sub-sample. Then, 
objects which are not detected above 3<r at 24 /rm/MIPS 
(i.e., Fv 24 < 20 — 25 pJy, on average) are indicated as MIPS- 
u galaxies, and the corresponding upper limits on their SFR 
are shown as red filled circles. Thus, this sub-sample should 
include quenched systems, or galaxies with a low residual of 
star-formation. The 17 secure AGN hosts identified in Sec¬ 
tion noi are highlighted by four-pointed stars, and the AGN 
candidate by a four-pointed star+circle. For comparison, we 
also show the parent sample of 24 ^m/MIPS-detected star¬ 
forming galaxies at 1.4 ^ z ^ 2, for which the IR contribu¬ 
tion to the SFR(IR+UV) was derived from Eq. ©), and (|5j 
(black open circles). 

As mentioned in Section im for completeness we also 
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included in our sample the objects with A 9k-irac > C//5, 
which are blended in IRAC and, consequently, in MIPS. A 
close examination of HST and MIPS maps for all the 56 
galaxies then revealed that 6 other objects are blended and 
could have been misclassified (i.e., as MIPS-d or MIPS-u) 
by the automatic procedure associating the MIPS sources to 
A-selected counterparts. Most of them reside in projected 
(or real) galaxy pairs, and in many cases both objects in a 
pair are included in our massive sample (see Appendix ®. 
Hence, we attempted to deblend these objects to derive a 
more accurate estimate of their SFR, and re-classify them 
as MIPS-detected or undetected, if necessary, as briefly sum¬ 
marized below and further detailed in Appendix IH). 

We used GALFIT to fit PSF models to the 24 pm 
maps with the F160W/WFC3 positions as priors, leaving 
the galaxy magnitude as the only free parameter to be recov¬ 
ered. This procedure worked quite smoothly in most cases 
but some “ambiguous” objects are discussed in Section Pf.d.ll 
In the right panel of Figure [2] we show the final SFR-M* 
diagram obtained after MIPS de-blending, while in the mid¬ 
dle panel we show the ‘migration path’ of each deblended 
galaxy from the position occupied before deblending (hor¬ 
izontal bar, color-coded as in the left panel) and after de¬ 
blending (filled circle, color coded as in the right panel). The 
‘migration path’ is also color-coded depending on the sub¬ 
sample to which the object belong after de-blending (e.g., 
migrations from the MIPS-d to the MIPS-u sub-sample are 
drawn in red, and vice-versa in green). From Figure[2]it ap¬ 
pears that several galaxies, previously classified as MIPS-d, 
have been reclassified as MIPS-u objects after deblending 
(MIPS-d—»u in Table [T]). This is the case for objects: #427, 
#880, #1187, #3230, #3231, and #5530. On the contrary, 
only galaxy #720, previously indicated as MIPS-u, has been 
reclassified as MIPS-d (MIPS-u—>-d). Then, object #7923 is 
confirmed to be a MIPS-u galaxy, but with a lower limit on 
its 24 pm flux, hence on its SFR. The de-blending procedure 
for each object is detailed in Appendix El 

All the PACS-d galaxies are within the MS band, 
while all the MIPS-u galaxies are more than 4x below it. 
The MIPS-d sub-class instead includes objects both on, 
and below the MS. Hence, some of the MIPS-d galax¬ 
ies could be quenching, or even quenched systems, judg¬ 
ing from the relatively wide range that they cover in sSFR 
(10 -9 — 10 -10 ' 4 yr _1 , i.e. from the typical MS rate down to 
~ 20x lower). It is also possible that some MIPS-d galaxy 
is instead undergoing a sort of rejuvenation process, caused 
by a merging/accretion event, although we expect that most 
such massive high-z galaxies are evolving in the opposite 
direction. Finally, in some case the observed MIPS emis¬ 
sion could be powered not (only) by star-formation, but 
(also) by obscured AGN , AGB stars, or old stellar popu¬ 
lation dSalim et al.ll2009h . 

In the following sections (cf. Section fl. 2114.31 and l4.3.Il) . 
we combined the IR/X/radio information with the results 
from SED fitting, and two-color diagrams to better constrain 
the SFR of our galaxies and shed some light on the nature 
of some ambiguous objects. 


4.2 SED fitting 

We built the multi-color SEDs by fitting the optical to 
IRAC/8 /im photometry using the hyperzmass software 


dBolzonella. Miralles fe Pellol |200C)I : IPozzetti et al.l l2007l l. 
For IRAC-blended objects we excluded the IRAC bands 
from the fit. However, we verified that the best-fit param¬ 
eters do n ot change substa ntially using instead the IRAC 
fluxes from I Guo et al.l (|2013l ^ . which are de-blended with the 
TFIT procedure. For galaxies at 1.4 ^ z ^ 1.7, the 8 gm 
flux was excluded from the SED fit, as at those redshifts it 
corresponds to a rest-frame part of the spectrum that may 
include thermal dust and 3.3 pm PAH emission line, features 
which are not included in the used stellar population models. 
We adopted BC03 stellar population models, with a grid of 
star formation histories (SFH) including an SSP, a constant 
SFR, and exponentially-declining r-models (r=0.1, 0.3, 1, 2, 
5, 10, 15, and 30 Gyr). To reduce the number of free param¬ 
eters, we only used templates with solar metallicity. These 
adopted star formation histories are easy to handle but are 
certainly inappropriate for near-MS star-forming galaxies at 
these redshifts, where actually the SFRs must have increased 
with time (on average) as a result of the mere existence of the 
MS jRenzin [2009|; IMaraston et al.ll 2010l: IPeng et a l.ll2010bl : 
iReddv et al. 20121 4 . As shown by IMaraston et al.l ( 201Ch . for 
star-forming galaxies, the age returned by the SED fit should 
be interpreted as the typical age of the stars contributing the 
bulk of the light, rather than the time since the beginning 
of star formation, as it formally represents in the fits. For 
this reason, to avoid estimating unrealistic low ages for such 
galaxies, in which the youngest stellar population outshines 
the older one, we imposed a lower age limit of 0.1 Gyr. 

The age returned by the SED fits has yet a different 
physical meaning in the case of quenched galaxies, i.e., here 
assumed to be those for which only upper limits can be 
placed on SFR(UV+IR). As apparent from Table [2] for 
many such galaxies the fits give r = 0, which is to say that 
a single age, simple stellar population (SSP) provides the 
best fit. Correspondingly, the age can be interpreted as rep¬ 
resentative of the time elapsed since quenching. This is a 
fair interpretation also for the other cases, in which r # 0 
but is << age, hence the current SFR. is very low in the best 
fit exponentially declining model. In the following, we adopt 
this interpretation for the ages of quenched galaxies. 

The SED fitting results are shown in Tableland Fig¬ 
ure [31 where symbols are like in Figure [21 In the left panel 
of Figure [3] we compare the best-fit stellar masses derived 
using BC03 models, and Salpeter IMF (M*sed,bco3), with 
those derived following D04 and D07, as explained in Sec¬ 
tion 13.11 (M* dim). The two quantities agree well, with an 
average ratio close to unity ((M*do4/A4+sed,bco3) = 1.05) 
and a Ict scatter of 0.16 dex. By inspecting the optical/nIR 
best-fit SEDs of the AGN hosts we found that most of 
them do not seems to suffer significant contamination by the 
AGN emission. On the contrary, two objects (i.e., #1906, 
and #6898) show a typical ‘IRAC power-law’ SED, indi¬ 
cating th at the AGN component is s wamping the stella r 
emission (lAlonso-Herrero et alJ 120061 : iDonlev et al.l 120071 ) . 
Hence, for these two objects, in particular, and also for all 
the AGN hosts with at least a MIPS detection (11/18), we 
verified that the adopted stellar masses (M*004) are consis¬ 
tent, (within ~ 0.25 dex) with those derived by fitting the 
optical to far-IR SEDs with a modified version of the MAG- 
PHYS softw are, which also includes the AGN warm dust 
comp onent (da Cunha, Chariot fe Elbazl 120081 : iBerta et al.l 
120131) . 
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Figure 2. In this figure we show the position of our massive sample galaxies on the MS at (z) = 1.7 (black solid line), before (left panel) 
and after (right panel) the MIPS-deblending procedure (see Appendix [7\l i . In the central panel we report the “migration path” of each 
deblended galaxy, from the position occupied before (horizontal bar, color coded as the object in the left panel), and after deblending 
(filled circle, color coded as the object in the right panel). The right panel shows the final positions with respect to the MS. The PACS-d, 
MIPS-d, and MIPS-u galaxies are displayed with different colors, as labeled in the figure, and explained in the text (Section [4j. Note 
that the red filled circles just represent 2, or 3 a upper limits to the SFR for MIPS-u galaxies (cf., Section l2.ll and 14.111 . For comparison, 
lower-mass star-forming galaxies at 1.4 ^ z ^ 2, and detected at 24 /vrn/MTPS. from the parent sample, are shown the as black open 
circles. 


While the stellar mass inferred from SED fitting is rel¬ 
atively robust, both for quen ched and star-forming galax¬ 
ies (cf. iMaraston et alj l20ld ). the other SED-fitting de¬ 
rived parameters are affected by the well-known degener¬ 
acy between the galaxy age and reddening (Av,sed), then 
also affecting the measured SFR(SED). However, if prop¬ 
erly combined with multi-wavelength data (i.e., IR, ra¬ 
dio, X), the SED fitting analysis can provide some in¬ 
sight on the nature of the studied galaxies. For instance, 
it can help to better constrain the reddening and SFR, 
for galaxies that are undetected in MIPS, and, in general, 
for galaxies below the MS. This is visible from the mid¬ 
dle and right panels of Figure [3] where the best-fit red¬ 
dening values and SFRs are compared with those derived 
from the IR and UV-uncorrected luminosities (Av,ir+uv 
and SFR(IR+UV)). For almost all the MIPS-u galaxies, the 
fit suggests lower values of Av,sed, and SFR(SED), with 
respect to the Av(IR-fUV) and SFR(IR-fUV) upper lim¬ 
its. This strengthen the interpretation of such galaxies as 
already quenched, or nearly so. When considering only the 
PACS-d and MIPS-d galaxies, Av.sed and Ay (IR+UV) are 
on average in good agreement, with a relative difference of 
(Ay ,sed — Av(IR+UV))/Ay(IR+UV))=-0.09 and lcr scat¬ 
ter of ~0.41. Some particular cases (flagged with numbers 
from 1 to 9) in the right and middle panels of Figure 0 are 
discussed in Section [4]3A] 


4.3 Observed and rest-frame two-color plots 

In Figure [4] we explore the observed and rest-frame colors of 
our sample in two color-color diagrams widely used in the 
literature to separate star-forming from passive galax ies at 
high redshifts, i.e., the BzK plot jDaddi et al.ll2004al ) and 
the rest-frame (U — V) vs. (V — J) (hereafter, UVJ) plot 


llLabbe et al.ll2005l : IWuvts et al.ll2007l : Williams et al.ll2009l ) . 

The different sub-samples are color-coded as labeled in the 
right panel, and detailed in Section 14.11 and Figure [7] The 
BzK diagram is a quite efficient criterion to separate galax¬ 
ies at 1.4 ^ z ^ 2.5 from those at lower-redshift and among 
these high -2 galaxies it separates the star-forming ones 
( sBzK , to be found in the left part of the diagram) from 
the passive ones ( pBzK , located within the wedge in the 
top-right side of the plot). Lower -2 systems are to be found 
in its bottom-right region of the plot. Rightwards arrows in 
the BzK diagram indicate lower limits to the B — z color, 
for objects undetected in R-band. These very red objects 
include both quenched galaxies (that would be classified 
as pBzK s with deeper R-ba nd data) and highly-red dened 
dusty star-forming galaxies (IRodighiero et al.1120141 ). The 
UV J diagram (right panel) has the advantage of allowing to 
separate dusty star-forming from quenched galaxies at any 
redshift. In fact, since the dust-free quenched objects have 
bluer V — J colors with respect to the extincted star-forming 
galaxies with similarly red U — V colors, they can be isolated 
on the top-left box of the diagram. We note that galaxies 
in this region define a tight sequence, which is almost par¬ 
allel to the diagonal border line of the quenched box. The 
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Figure 3. The SED fitting derived Mass, SFR, and reddening are compared with the same quantities, as derived with different methods. 
In all the panels symbols are like in Figure [2] The numbers from 1 to 9, labeled in the middle and right panels, represent the odd cases 
discussed in Section 14.3.11 and correspond to objects: #428=1, #557=2, #720=3, #856=4, #880=5, #2940=6, #3853=7, #5509=8, 
#6193=9. Left panel: Best-fit masses (M*sed) us masses derived from ff-band magnitude (Mwxm). Middle panel: Best-fit reddening 
(Av.SBd) us A v (IR+UV), derived from the ratio between the total and UV uncorrected SFR. Right panel: Best-fit SFR (SFRssp)us 
SFR(IR+UV), derived as explained in Section 13.21 


extent of the “quenched sequence” mostly depends on the 
age range spanned by this population, which in our sample is 
~ 0.5 — 2 Gyr. In fact, young and old quenched objects could 
be split according to the ir V — J color along the sequence 
(cf. iWhitaker et a.1.1 l201Clll . The tightness of the “quenched 
sequence” in the UV J plane is instead related to the fact 
that most of these systems are affected by a small amount 
of dust reddening. The U — V color of star-forming galaxies 
depends on their dus t extinction, and sepa rates blue from 
red dusty objects (cf. iFumagalli et al.ll2014l 'l. 

In the left panel of Figure[X]some objects do not match 
the BzK criterion, i.e., occupy the lower part of the plot 
that should contain z < 1.4 galaxies (namely, #282, #427, 
#686, #887, #1187, #1272, #3230, #4705, #5556, #6193, 
#7491, #7617, and #7923). For some objects (i.e., those 
closer to the boundaries of the pBzK or sBzK regions) this 
could be due to even small photometric errors, and in one 
case to Type-1 AGN contamination (#686). Then, 4 of the 
7 MIPS-u galaxies with the bluest z — K colors (i.e., #282, 
#7491, #7617, and #7923), do not have spectroscopic iden¬ 
tification, and could be at slightly lower r edshift. However, as 
shown in Figure 8 of lDaddi et al.l (l2004blh the colors of these 
objects are also consistent with those of relatively young 
(< 2 Gyr) passively evolving galaxies at z ~ 1.4 — 1.5, char¬ 
acterized by a rapidly declining SFR, which could be missed 
by a strict BzK selection. It is worth noticing that most of 
these bluer MIPS-u galaxies also lie in the young-quenched 
region of the UV J diagram, and have SEDs best-fitted with 
SSPs, or declining r = 0.1 Gyr models with age ~ 0.7 — 1.3 
Gyr. Hence, we did not exclude any of these sources from 
the sample. 

In Figure [4] the galaxy colors strongly correlate with 
their position relative to the MS: MIPS-u galaxies are gen¬ 
erally in the “quenched regions” ( pBzK , and pUVJ), and 
PACS-d galaxies in the star-forming regions ( sBzK , and 
sUVJ), in both the panels. With two exceptions, the MIPS- 
d objects lie in the star-forming region in the BzK diagram, 
while half of them are in the quenched region in the UV J 


diagram. In the next section, we will discuss these objects, 
together with other peculiar cases. 

4-3.1 A few odd cases 

A few objects show BzK or UV J colors and best-fit SED 
which are apparently at odd with respect to their position 
relative to the MS. In some cases this is due to (even small) 
photometric errors, but in other cases this could indicate a 
different nature of these “outliers” (i.e., star-forming galax¬ 
ies with a relatively old age, on the way to be quenched, 
or rejuvenated, AGN, etc.), as discussed below. To facili¬ 
tate the identification of these objects in Figure [3] and [4] we 
flagged them with numbers from 1 to 9. 

Object #428 (‘1’ in the figures) is a MIPS-u galaxy, 
but shows colors similar to star-forming galaxies, i.e., lies in 
the sBzK and sUV J regions. Judging from its appearance 
(see Figure [9J, this object may be a lens Early-Type galaxy 
surrounded by a blue arc. So, the blue colors can be due 
to the magnified background source when the photometry 
does not separate the arc from the central galaxy. Indeed, the 
rest-frame colors of the central galaxy, having masked out 
the blue arc as detailed in Section fo~2l are (U - V) = 1.93 
and (V — J) = 1.58, consistent with those of a quenched 
object (see Section liPl) . 

There are two MIPS-u objects in the dusty-star-forming 
box of the UV J plane. They are #880, and #5509 (flagged 
as ‘5’, and ‘8’), and in the BzK diagram are shown as lower 
limit to B — z, being undetected in R-band. Their BzK 
colors are compatible with those of old red galaxies, but 
the SED fitting analysis indicates a higher reddening and 
SFR(SED) with respect to other MIPS-u objects (cf. Fig- 
ure[3]and Table[2j, yet still in agreement with the upper lim¬ 
its derived from the 24 /rm/MIPS. In particular, since object 
#880 is one of the MIPS-deblended sources (see Section l2Tl 
and Appendix E} , we do not exclude that it has a low resid¬ 
ual star-formation (partly contributing to the MIPS blended 
source). However, it is also very close to the boundaries of 
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Figure 4. BzK and UVJ diagrams for the three sub-samples, where MIPS-u are the red filled circles, MIPS-d are green filled circles and 
PACS-d are blue filled circles, as labeled in the right panel of the figure. X-ray sources and IRAC power-law are also highlighted by cyan 
stars. The odd cases discussed in Section 14.3.11 are flagged with numbers from 1 to 9, as in Figure[3] 


the pUV J region, and, alternatively, its odd position in the 
UV J color diagram could just depend on small photometric 
errors. The AGN host #5509 is undetected, or very faint, in 
all the optical bands. Its location in the dusty-sf/FJ could 
be due to residual obscured star-formation, which could also 
be responsible of the marginal 24 /rm/MIPS detection of the 
object (2.5cr, see Section [4] and Table [I}. It is also possible 
that the AGN partly, or totally, contributes to the low MIPS 
emission of the source, and that its odd position on the UV J 
diagram is instead due to a hig her redshift w ith respect to 
the photo-,z qu oted here dGrazian et alj|200d . for instance, 
ILuo et ai1l2010l and lHsu et al.ll2014l report z p hot — 2.5). 

Object #2940 (flagged as ‘6’) is a MIPS-u AGN host 
lying in the sBzK region. It is in the quenched region of 
the UVJ diagram, but its bluer B — z color could indicate 
that a low residual star-formation is still present. The SFR 
inferred from the best-fit SED is ~ 8 Mq/ yr, which is very 
close to the upper limit from MIPS (cf., the right panel of 
Figure [3j. 

Two of the MIPS-d galaxies > 4x below the MS, i.e., 
#557 and #856 (labeled as ‘2’, and ‘4’, respectively) lie 
in the quenched region, both in the BzK , and in the UVJ 
diagram. Moreover, their best-fit SEDs indicate a very low 
SFR (i.e., SFR(SED)<3 M 0 /yr), and modest extinction (cf. 
Tableland Figure^. Hence, they are likely to be quenched 
galaxies, where the 24 /rm/MIPS detection (slightly above 
3a) could be due to old AGB stars, or alternatively to a 
highly obscured AGN, rather than to star-formation. In the 
following we consider them as quenched galaxies. 

On the other hand, the third MIPS-d galaxy below 
the MS, #3853 (labeled as ‘7’), seems to have a differ¬ 
ent nature with respect to the other two discussed above. 
It also lies in the quenched region of the UVJ diagram, 
and has a red, z — K color (~ 2.4), consistent with a rel¬ 
atively old age (1.14 Gyr, cf. Table [2j, but it is classified 


as a sBzK, due to the bluer B — z color. This seems to 
indicate that the 24 /rm/MIPS emission is plausibly pow¬ 
ered by a low residual star-formation, as also confirmed by 
the SED fitting results in Figure [3] {Ay (SED) = 0.99, and 
SFR(SED)=20 Mg/yr, in agreement with that derived from 
MIPS, SFR(IR+UV)=23±6M 0 /yr ). 


The MIPS-d AGN objects #720 and #6193 (i.e., ‘3’ 
and ‘9’) have observed colors typical of sBzK but lie in the 
“young-quenched” region of the UVJ diagram. In fact, they 
are redder with respect to the other MIPS-d galaxies on 
the MS, and best-fit with declining SFR with r = 0.1 stel¬ 
lar population models, age=0.64 Gyr, low reddening, and 
SFR(SED)~ 2 — 3 M 0 /yr. These results do not agree with 
the high SFR(IR+UV)> 100 Mq/ yr and large reddening 
(~ 1.8 — 2 mag) derived from the 24 /rm/MIPS emission of 
these galaxies. If the AGNs were responsible of most of the 
MIPS emission, both of these objects should be shifted be¬ 
low the MS in Figure[2] entering the sub-sample of quenched 
galaxies. However, it is not trivial to estimate the fraction 
of the IR luminosity due to the AGN and to obscured star- 
formation, because both objects are undetected in the far- 
IR. The bluer observed colors, and the younger best-fit age 
of these galaxies, with respect to other MIPS-u quenched 
galaxies in the sample (including AGN hosts), seems to sug¬ 
gest that they are still star-forming, but maybe on the way 
to be quenched (possibly triggered by the AGN). On the 
contrary, the third AGN among the MIPS-d galaxies, i.e., 
#552, shows BzK/UVJ colors and best-fit SED consis¬ 
tent with a moderate unobscured star-formation activity, in 
agreement with what expected based on the MIPS emission 
(SFRsbo r;SFR(IR+UV)« 56 M 0 /yr). 
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4-3.2 Summary of the source classification 

As shown in the previous sections, for galaxies lying more 
than 4x below the MS the best-fit SEDs provide stricter 
upper-limits to their SFR, compared to the limits set by the 
24 /rm/MIPS data. Hence, hereafter we use these new upper- 
limits to the SFR for all the MIPS-u galaxies, as well as for 
the three MIPS-d galaxies located below the MS (i.e., #557, 
#856, and #3853, flagged as ‘2’, ‘4’, and ‘7’ in Figures [3][I]). 

Figure [5] shows the distribution of the sources as a func¬ 
tion of the “distance” from the main sequence, parametrized 
as the sSFR normalized to that predicted by Eq. JT} for 
a MS galaxy of the same mass, in logarithmic scale, i.e. 
log(sSFR/sSFRivis)). The red solid and dot-dashed lines in 
the figure show the average MS, and the offsets of ±0.6 dex, 
corresponding to 4x above, and below it. The MIPS-d galax¬ 
ies (undetected in the FIR, hence which are not PACS-d 
galaxies) are highlighted with a white frame. Within the MS 
band all other galaxies are detected in the FIR, while outside 
this band none of the galaxies is detected. AGN hosts are 
identified by a cyan frame. The first log(sSFR/sSFRMs) bin 
of the histogram includes only galaxies which are best-fit by 
an SSP, hence with SFR(SED)=0. For each of the remaining 
bins, the postage stamps are displayed in such a way that the 
distance from the MS also increases, i.e., log(sSFR/sSFRMs) 
decreases, from top to bottom along the y axis. 

Out of the 56 most massive galaxies at 1.4 ^ z ^ 2 
in GOODS-S, 25 objects (~ 45%) are on the MS, i.e., 
within a factor of 4 (±0.6 dex in Figure [Sj from the av¬ 
erage MS relation at z — 1.7, while 31 objects (~ 55%) 
are below it (log(sSFR/sSFRMs) <-0.6). Several authors 
have found a flatteni ng of the Main Se q uence slope at the 


high-mass end (e.g., 

Bauer et al.l 1201 ll: IWuvts et al.l 1201 ll: 

Whitaker et al. 20121: 

Rodiehiero et al. 2014), interpreted as 


a further evidence that the mass-quenc hing process is mor e 
efficient on the most massive galaxies ( IPeng et al.M2010bi l. 
Also here the fraction of MS galaxies below the average re¬ 
lation (2/3) exceeds the fraction above it (1/3), although the 
sample is too small to establish the significance of this re¬ 
sult. On the other hand, if the transition from star-forming 
to passive galaxies is linked to a morphological transfor¬ 
mation, one cannot exclude that the appearance of passive 
bulges in the most evolved massive star-forming galaxies, 
contributing to the galaxy M* but not much to its SFR 
could flatten the MS slop e at the highest masses (lLang et al.l 
120141 ; iBruce et al.ll2014 and Tacchella et al. 2014, submit¬ 
ted). Actually, a visual inspection of the image stamps in 
Figure 0 seems to confirm that many MS galaxies with 
—0.6 <log(sSFR/sSFRMs) < 0 host a bulge. As shown in 
the following, this has been also confirmed through a 2D 
surface brightness fitting in the HST/ WFC3 images (Sec¬ 
tion o. 

At the end of this exercise, i.e., after the de-blending 
procedure, and the use of SED fitting to better constrain 
the SFR for galaxies below the MS, we conclude that all the 
galaxies below the MS band have a sSFR more than 10 x 
lower than that of MS galaxies in the same redshift range, 
and most of them ~ 100 to ~ 1000 times lower (Figure [5]). 
Their global properties are also typical of quenched objects. 
This suggests that at z ~ 1.4 — 2 the mass-quenching may be 
a relatively fast process, much shorter than the time spent 
by these galaxies in their massive phase within the MS band 


(i.e., a transition time <<~lGyr), as one may expect if 
quenching is caused by the sudden expulsion of all the gas 
and dust from the galaxy. 

In summary, our sample of 56 galaxies splits in two 
main groups: (i) the MS galaxies, included in the MS band 
(—0.6 <log(sSFR/sSFRMs) < 0), and (ii) the quenched 
galaxies, with log(sSFR/sSFRMs) < — 1, which have a 
mass doubling time (the time needed by the galaxy to 
double its M*, in case it continues to form stars at the 
present rate) ^ 10 Gyr. In our small sample a well-defined, 
sizable class of objects on the way to be quenched can¬ 
not be unambiguously recognized. Such transition galax¬ 
ies could be found among the MS galaxies with the low¬ 
est sSFR. (—0.6 <log(sSFR/sSFRMs) < 0), and/or among 
the quenched galaxies with the highest sSFR upper-limits 
(—2 < log(sSFR/sSFRMs) < —0.6). Clearly a much bigger 
sample should be analyzed to adequately populate this class 
of objects which must exist and should allow to estimate 
the quenching timescale. However, the present investigation 
demonstrates that a careful, object-by-object study is nec¬ 
essary to identify bona-fide transition objects, given that 
blending may lead to gross errors in the estimated SFRs, 
especially for this class of galaxies (see Figure [2j. 


5 MORPHOLOGICAL CLASSIFICATION 
5.1 Single-Sersic surface brightness fitting 


We used GALFIT to fit the surface brightness (SB) profile 
of our galaxies. Since the sample includes both star-forming 
and quenched high-z massive galaxies, for a first analysis we 
fit the SB profiles of each galaxy in the F 160W/ WFC3 filter 
(77-band) with a single-Sersic function dSersi d 1196811 . leav¬ 
ing the Sersic index (n) as a free parameter. I11 general, the 
surface brightness distribution of high-z galaxies is far from 
following a regular Sersic profile, and especially so those of 


star-form ing galaxies, that often are clumpy or merger-like 
systems (|F5rste r Schreiber et al] 120091 : lEpinat et all 120091 : 
iGenzel et al. 20111 1. However, this method allows us to clas¬ 
sify objects as either Early-Type Galaxies (ETG, bulge- 
dominated, n > 2 — 2.5) or Late-Type Galaxies (LTG, disk- 
dominated, n < 2 — 2.5). Three out of the 56 galaxies have 
irregular morphology, and show strong residuals when fit¬ 
ted with a single Sersic profile, or with a bulge±disk model. 
They are the AGN ‘candidate’ #5503 (see Section [373]) , and 
the inactive galaxies #3236, and #5534, which have been 
excluded from the morphological analysis presented in the 
following. To reliably measure the sky background, we run 
GALFIT on large images, of 30 arcsecond in size centered 
on each galaxy, and fit together all the brightest sources 
(H <23.5) included in the field ( masking out the fainter 
ones). A prior run with SExtractor dBertin fc Arnoutsll 19961 ) 
on the same images provided the morphological parameters 
to be set as starting point in the fit. The PSF model was 
built using the 10-20 closest unsaturated stars in the field. 
All the parameters, i.e., sky background, Sersic index (n), 
centroid, total magnitude (nii6o), half-light semi-major axis 
(se), position angle (PA), and axial ratio (6/a), were left as 
free. The results are reported in Tableland Figure [6] 

The left panel of Figure [6] shows the index n as a 
function of the distance from the MS at (z) = 1.7, i.e., 
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Figure 5. The distribution of the sample galaxies as a function of the ‘distance’ from the main sequence, i.e., the log(sSFR/sSFRjv[s))- 
As explained in the text, hereafter, we used the SFR inferred from the best-fit SED (SFR(SED)) for all the galaxies below the MS strip. 
Within each logfsSFR/ sSFRms )) bin the distance from the MS also increases from top to bottom. On the contrary, in the first bin, which 
only includes galaxies best-fit with a SSP, and with SFR(SED)=0, objects are ordered randomly. The red solid and dot-dashed lines 
show the MS, and the loci 4x above and below it, respectively. The MIPS-detected (and PACS-undetected) galaxies are distinguished 
by the white frames, and AGN hosts by a cyan frame. The postage stamps are red-green-blue (R-G-B) composite images (R=F160W, 
G=F850LP, B=F435W) of 3 arcsec size. 


log(sSFR/sSFR M s)• The PACS-d, MIPS-d, and MIPS-u 
sub-samples are shown in blue, green, and red, respectively, 
and the AGN highlighted with cyan four-point stars, such as 
in Figures [2141 The vertical solid, and dash-dotted lines rep¬ 
resent the MS (Eq. @), and its 4x lower boundary (labeled 
as “MS/4”), respectively. The Sersic n = 2.5 line shows the 
separation between “Early-Type” (or bulge-dominated) and 
“Late-Type” (or disk-dominated) galaxies. The right panel 
of the same figure, shows the distance from the MS vs. the 
distance from the relation of local ETGs (i.e., the 

circularized effec tive radius, R e , n ormalized to the value at 
z — 0, Ro, from lShen et al.ll2003l 'l. Here we use the M,-R e 
relation of local ETGs as a reference for both the quenched 
and the star-forming galaxies in our sample, as z ~ 1.4 — 2 
star-forming galaxies with such high masses (and SFRs) 
will necessarily be quenched in a relatively short time (to 
avoid mass overgrowing), and then passively evolve down 
to 2 = 0 without major further increase in mass. The bot¬ 
tom panels show the mean values of the same quantities, in 


log(sSFR/sSFRMs) bins chosen so to have 7—11 objects per 
bin. Black filled circles are relative to all the objects within a 
bin, while red open squares are obtained by only considering 
galaxies without evidence of AGN. 

From Figure [6] it appears that almost all the quenched 
galaxies show an ‘early-type’ morphology, with n > 2 — 2.5, 
and a circularized effective radius ( R e ) from 1.4 to 5 (2.5 
on average) times smaller with respect to the local ETGs of 
the same mass, a well known property of high -2 quenched 
galaxies. On the contrary, MS galaxies split in two groups, 
above and below the border-line at n = 2.5. Almost all the 
MS galaxies with n > 2 — 2.5 host an AGN, whereas all the 
AGN host have n > 2 (but the AGN ‘candidate’, #5503, 
which has irregular morphology and is not shown in the fig¬ 
ure). The high Sersic index of the AGN hosts could be due 
either to the nuclear emission of the AGN itself, detected as 
a central, unres olved point-source, or to a significant bulge 
component (cf.. IRosario et aDl2015l '). We investigate this is¬ 
sue in Section 15.21 by means of surface brightness profile 
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log(sSFR/sSFR|^) log(sSFR/sSFR^g) 


Figure 6. Left: Sersic index vs distance from MS (i.e.,log(sSFR/sSFRjvfS'))- The vertical red dot-dashed line at log(sSFR/sSFRM 5 ) = 
—0.6 marks the lower boundary of the MS band, and the horizontal n = 2.5 line separates disk-dominated and bulge-dominated objects. 
Right: The “distance” from the M t -R e relation for local ETGs (horizontal solid line), expressed as R e /Ro, as a function of the distance 
from MS, log ( s S F R,/ s S F Rm S )■ Symbols are like in the previous figures. Bottom panels: Mean values of the same quantities, in 
Iog(sSFR/sSFRvig) bins chosen so to have 7 — 11 objects per bin. Black filled circles are relative to all the objects within a bin, while 
red open squares are obtained by only considering galaxies without evidence of AGN. 


decomposition, in bulge+disk, and bulge+disk+PSF com¬ 
ponents. 

MS galaxies have in general size comparable with those 
of local ETGs, a part from 4 AGN hosts (i.e., #1906, #6898, 
#6193, #6352), and one normal galaxy (#5415), which are 
comparable t o the compa c t star - formi ng galaxies (CSFGs) 
identified by iBarro et al. ( 120131 . 120141 . hereafter, B13 and 
B14) at 1.4 < z < 3, and considered the direct progeni¬ 
tors of compact quenched systems at similar redshifts. Since, 
as mentioned in Section 14.21 #1906, and #6898 are con¬ 
taminated by the AGN in their photometric SEDs, it is 
not clear if their compact size are intrinsical, or just due 
to the AGN affecting the surface brightness profile (cf., 
also the discussion in Section [Oil . In any case, all these 
5 objects have a relatively low sSFR., lying in the region 
—0.6 <log(sSFR/sSFR,Ms) < 0 of Figure [6] consistently 
with having started their quenching phase. 


Concerning the quenched population, our results are 
in agreement with the well-documented smaller size of 
the high -2 quenched galaxies with respect to th e local 


ones 


(e.g. Daddi et al. 20051; Truiillo et al.l 20061. 

20071; 

Longhetti et al. 2007);!McGrath et al. 20081: 

Toft et al. 

20071 

Cimatti et al. 20081: Ivan Dokkum et alJ2008 

: Buitrago et al. 

2008 


Damianov et al. 20091: Saracco. Longhetti & Andreon 

200S 

Cassata et al. 20101: Barro et al. 2013: Bezanson et al. 

2013 

: 

Szomoru et al.l 20131). On the other hand, not all the 


quenched galaxies at high -2 have been foun d to be under¬ 
sized, both in the field and in proto-clusters ( Mancinij^_aL 


201C; 

Saracco. Longhetti & Gargiulo 201 ll: Onodera et al. 

2012; 

Stott et al. 2011: Huertas-ComDanv et al. 20131). It 


has been suggested that the average size growth of 
ETGs from g = 3 to g = 0 is not (or not only) 
due to the evolution of the individual galaxies (e.g., by 


mergers), but mostly (or partly) caused by the grad¬ 
ual additi on of larger, newly - quen ch ed galaxies at lower 
redshifts ([Vale ntinuzzi et alf l20ld : ICassata et al.l 120131 : 
IPoggianti et al. 20131 : Carollo et al.ll2013l). although some 
authors have a different view (e.g., van der Wei et al.ll2009l ). 
In our sample, the scarcity of compact with respect to nor¬ 
mal star-forming galaxies, could be interpreted as an indi¬ 
cation that the majority of the incoming quenched galaxies 
at 2 < 1.4 are progressively larger in size with respect to the 
objects quenched at higher redshifts. However, the small size 
of the sample prevent us from providing new clues on this 


5.2 Bulge-Disk decomposition 

We built two-component bulge+disk models, by fitting each 
galaxy with two Sersic components with n = 1 (disk) and 
n = 4 (bulge), both centered on the position derived in 
the previous single-Sersic fit. The bulge to total flux ratio 
(hereafter, B/T) gives an indication of the relevance of the 
bulge-component in the galaxy. We discarded those mod¬ 
els in which unreliable parameter values were retrieved for 
some of the subcomponents, i.e., too large/too small effec¬ 
tive radius (R e > 4"0 or R e < OVOS), or axis ratio b/a < 0.1. 
In most cases, the component with unreliable parameters 
was als o found to contrib ute less than 10% to the total 
flux (see lBruce et al.ll2012l . for a similar discussion). In sum¬ 
mary, beyond the three objects with irregular morphology, 
the double-component model was discarded for 15 galaxies, 
of which 10 are classified as simple disks (B/T=0), and 5 as 
simple bulges (B/T=l), (see Tabic [3]). In Figure[7]we show 
the relation between the B/T and the single-Sersic n param¬ 
eter for our sample (color-coded as in the previous figures). 
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Since [Bruce et ail (I2OI2I ) also studied bulge+disk decompo¬ 
sition for similarly massive galaxies (M* > 1.7 x 10 11 Mq, 
Salpeter IMF) at 1 < z < 3, for comparison, we also show 
in the figure objects included in the same redshift and mass 
ra nge from that sample (black small points). As in the 
iBruce et alj (l2012[ j sample, also for our galaxies there is a 
good correspondence between the morphological classifica¬ 
tion based on the single-Sersic index, n, and on B/T, al¬ 
though with a considerable scatter. In general, objects with 
B/T ^ 0.5 show n ^ 2 — 2.5, and those with n < 2 — 2.5 
show B/T < 0.5. However, there are also some objects with 
a high n and B/T< 0.5, i.e., located in the bottom-right 
quadrant of Figure^ Some of these objects have a relatively 
high 0.3 < B/T < 0.5, still consistent with the presence of 
a significant bulge, as suggested by their high Sersic index. 
For the remaining objects (two PACS-d galaxies, i.e.,#3066, 
and #6572, and three MIPS-u galaxies, i.e., #1084, #1187, 
#8121) the apparent incongruity between n and B/T may 
be due to the presence of an unresolved point-like source in 
the center, interpreted as a small bulge in the 2-component 
fit (.ReBuige < 1 pixel), that could also be due to residual 
star formation or an AGN in the center. This alternative 
is supported by the fact that when the surface brightness 
profiles of these galaxies are fit with a single-Sersic model 
+ central PSF, the resulting Sersic index is n ~ 1 (a pure 
disk), and the point-like component is found to contribute 
> 10% of the total galaxy flux (see Tabled]). Then, #428 
(flagged in the figure) is the already mentioned lens galaxy 
candidate (also highlighted in Figure[4] see also Section r5.3l) . 
In this case the low B/T could be explained by the presence 
of the magnified background source, in the outer part of the 
galaxy, fit as a relatively bright disk. 

The top panel of Figure [5] shows the bulge to total light 
ratio B/T vs log(sSFR/sSFR.Ms) in logarithmic scale, and 
the bottom panel the average trend, {B/T), as in Figure [6] 
As mentioned for the n index, most of the quenched galax¬ 
ies have B/T > 0.5, while the MS galaxies again split in 
two groups with B/T < 0.5 (64%) or > 0.5 (36%) (when 
considering a boundary of B/T = 0.3 they split almost in 
half). For galaxies within the MS band {B/T) seems to in¬ 
crease with decreasing log(sSFR/sSFRMs) (from ~ 0 — 0.2 at 
log(sSFR/sSFRMS > 0, to 0.5-0.7 at log(sSFR/sSFRMS < 
0), and then to remain almost constant below it, where al¬ 
most all the galaxies are bulge-dominated. Figures [7] and [8] 
show that the high single-Sersic index measured for all the 
secure AGNs on the MS appears to be due to the presence of 
a significant bulge component {B/T ^ 0.3—0.5) in almost all 
of them (9/10). The only AGN with B/T < 0.3 is #3066, al¬ 
ready discussed above. Since for most of the X-ray detected 
AGNs in our sample the nuclear emission does not dominate 
the photometric light in the optical-near-IR SED of the host 
galaxy, we expect a similarly low contamination also in the 
surface brightness profile. However, to rule out the possi¬ 
bility that the nuclear emission have biased the results to¬ 
wards higher Sersic indices and B/T ratios, we also tried to 
fit each object by adding a central point-source to the dou¬ 
ble component fit (i.e., bulge+disk+PSF). For most of the 
galaxies (50/56) the hypothetical point-source was found to 
contribute < 10% to the total flux. Moreover, the B/T ratio 
derived from the three-components fit (bulge+disk+PSF) 
agrees very well with that derived from the two-components 
fit (bulge+disk), for all the objects (including the X-ray 



Figure 7. Left: Bulge to Total light ratio (B/T) vs n Sersic . 
Symbols are like in the previous figures. The IBruce et all (i 201 2i ) 
results for objects in the same mass and redshift range of our 
sample are also shown for comparison (small black dots). 


detected AGNs). This means that the nuclear emission, if 
present, does not appreciably affect the light profile of our 
galaxies. The only exception may be represented by the two 
TRAC power-law’ sources, i.e., #1906, and #6898, already 
discussed in Section na and EH Their compact size, high 
Sersic n, and the fact that unreliable parameters were re¬ 
covered for the disk component, both in the Bulge+Disk, 
and in the Bulge+Disk+PSF fits (B/T=l), seems to sug¬ 
gest that their light profile is outshined by the nuclear emis¬ 
sion. Hence we cannot reliably estimate their morphological 
parameters. However, for all the other AGNs we conclude 
that a real link should exist between the presence of an 
AGN and a relevant bulge component in the host galaxy. 
This fits well with the notion that galaxies and supermas- 
sive black holes (SMBH) grow together in order to estab¬ 
lish the SMBH- bulge mass relation as observed in the local 
Universe (e.g.. [ Magorri an et al.l Il998l : iFerrarese fe Merrittl 
120001 : [Daddi et al.l 2007a ). The above analysis confirms that 


MS galaxies with M* ^ 10 11 Mq show on average a 
more relevant bulge component with r es pect to their lower 
mass counterparts (lLang et al.l 120141 : [Bruce et al.l 120141 : 
iTacchella et al.l 120141 ) . Moreover, since the most massive 
among star-forming galaxies are the first to be quenched, 
this supports the idea that the mass quenching processes 
is strongly related to the bulge growth and to AGN activ¬ 
ity, as a lso assumed by some theoretical and semi-analy t ical 
mode ls (ISilk fe Reedl 19981: iDi Matteo. Sprineel fe Hernauistl 


120051 : ISomerville et al.' 120081 ). This scenario is further sup- 
ported by the large fraction of AGNs that we found also 
among quenched galaxies, as discussed in Section [6] 
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Figure 8. Bulge to Total light ratio vs distance from MS, 
log(sSFR/sSFR]y[s)- Symbols are like in the previous figures. 
Bottom panel: Mean B/T ratio in bins of log(sSFR/sSFR]yis)) 
as in Figure [6] Black filled circles are relative to all the objects 
within a bin, and red open squares are obtained by excluding the 
AGNs. 


5.3 A double lens system? 

The close pair of galaxies #428 and #427 represents a puz¬ 
zling system as illustrated in Figure [9] and we suspect that 
they may form a double lens system. The red, compact 
galaxy #428 is surrounded by a blue arc, which is highly 
reminiscent of a section of an Einstein ring, as illustrated in 
the left panel of the figure. If this is indeed a lensing sys¬ 
tem, the very massive ETG #427, may contribute to further 
distort the image of a background source. 

Object #427 at z spe c = 1.427 is spectroscopically clas¬ 
sified as a passively evolving g alaxy based on the C aH+K, 
Mgll, and Mgl absorption lines (IVanzella et al.|[2008l . see Ta¬ 
ble 1), while no spectroscopic confirmation exists for #428 
which has a z p hot= 1.59. The small projected distance (1/88, 
i.e., ~16 kpc at this redshift) and the massive nature of the 
two galaxies makes it possible that they are a gravitation¬ 
ally bound pair at z = 1.427. The central spheroid of #428 
is undetected in the 17-band and close to the detection limit 
in the HST /ACS B-band, while the outer arc is detected in 
all bands, hence it is not a 17-band dropout. If indeed lensed, 
the source should be at 2 < z < 3. The radius of the blue 
arc is Bare = O'/S, somewhat larger than expected for the 
Einstein radius of an isolated Mtot ~ 40 x 10 11 Mq lens at 
z = 1.427 and a z ~ 2 — 3 background source, for which 
the Einstein radius would be Re = 0/33 — 0/6. This is as¬ 
suming a dark matter contribution ~ 40 times the stellar 
mass of the galaxy. Hence, the unusual configuration of the 
system, with only half visible ring around object #428 (and 
such large Re), could be due to the presence of the mas¬ 


sive #427 neighbor perturbing the gravitational field and 
producing a warped and atypical source image. 

The further exploration of the various components of 
this system has been pursued by fitting the surface bright¬ 
ness distribution of #428 and #427 with both single-Sersic 
function and with bulge+disk models (Section [SJ. From the 
single-Sersic fit we obtained a high Sersic-index (n = 3.8) 
and a relatively large Retire (5.8 kpc) for #428, and n = 
2.73 and Retire = 2.4 kpc for #427. The corresponding 
model is shown in the middle panel of Figure [jj] Once the 
single-Sersic models have been subtracted, strong ring-like 
residuals are visible in the outskirts of both galaxies as 
shown in the right panel of Figure [5] From this figure it 
also appears that the bottom-right part of the blue ring 
of #428 emerges in the residual image, with an evident 
brighter “blob” SE of the galaxy. A smaller ring also ap¬ 
pears in the residual image of #427, along with a struc¬ 
ture apparently bridging these two rings. The bulge+disk 
decomposition resulted in a very faint and small bulge 
( B/T ~ 0.2, Retire,bulge ~ 0.5 kpc) for #428, and a bright 
and large disk R e ,circ,disk ~5.2 kpc), but produced similarly 
strong residuals. This shows that an exponential disk is not 
suitable to fit the outer structure in #428. 

Thus, from these evidences we cannot conclude whether 
we are dealing with a double-lens system, or whether this 
complex structures are the result of tidal interactions. What 
is more certain is that the blue ring is a distinct, star-forming 
galaxy, either being lensed or torn apart by #428. We then 
derived more reliable morphological parameters and total 
photometry for the supposed lens galaxy by masking-out on 
the HST 77-band image the outer arc with an annulus of in¬ 
ternal and external radii of n n = 0/6 and r ou t = 1/2, respec¬ 
tively. By fitting again the galaxy with a single-Sersic func¬ 
tion we then derived nearly the same Sersic index, (n ~ 3.8), 
but a smaller effective radius, 7? e ,circ=3.9 kpc. Using the 
same masking we have also determined the total magnitudes 
in all the other available HST optical (i.e., B,V,i,z) and 
near-IR. ( Y,J ) filters and all the other GALFIT parameters 
to those derived in the 77-band. As shown in the Appendix 
(Figure |ATJ, a blended detection at MIPS/24 /im is also as¬ 
sociated to the system with the MIPS peak being clearly 
shifted toward S-E with respect of both the IRAC and the 
HST/WFC3 centroids. Hence, the 24 pm flux most likely 
comes from the blue arc, rather then from either of the two 
galaxies #427 and #428, confirming their quenched nature. 

We fit the new HST optical/near-IR photometry with 
an SED including relatively low reddening (Ay < 1), as 
inferred from the MIPS upper limit (assuming MIPS non 
detection, cf. Section 14.21 and Figure [3]). The derived sSFR 
(log(sSFRfsSFRMs) = —1.4) is consistent with that of 
a purely quenched galaxy. Whatever the nature of the 
#427/#428 pair, either a lensing or a tidally interacting 
system, it clearly deserves further study. 


6 AGN AND QUENCHING 

Out of the 56 massive galaxies in this study, 18 (or 17) 
are hosting an AGN, as revealed by their X-ray and ra¬ 
dio detection (cf. Section 3.3). This fairly high ~ 30% frac¬ 
tion is not unc ommon among massive galaxies at th ese high 
redshifts (e.g.. lDaddi et al.ll2007al ; iFiore et al.ll2009l j. AGN- 
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Figure 9. From left to right we show : the RGB image (where R=H-band (F160W), G= 2 -band (F850LP), B-band (F435W) of HST), 
the GALFIT model in the F160W filter for single-Sersic fit, and the residual (image-model) in the same filter for objects #427 and #428. 


driven, high velocity nuclear outflows are found to be ubiq¬ 
uitous among massive (> 10 11 Mq) star-forming, near-MS 
galaxies (Forster Schreiber et al. 2014; Genzel et al. 2014). 
Given their extension (few kpc) such outflows must have 
fairly long lifetimes (> 10 8 yrs) and therefore be detectable 
even if the powering AGN is caught in an off phase, due to 
its short term variability. This may provide further evidence 
for SMBHs being responsible for regulating the growth of 
the galaxies and vice versa. It also provides circumstantial 
support to the notion of AGN feedback as a plausible phys¬ 
ical process leading to mass quenching. However, evidence 
of AGN feedback, e.g., in the form of fast nuclear outflow, 
is not necessarily evidence for AGN quenching, as pointed 
out in the introduction. 


Particularly interesting is the fact that the AGN frac¬ 
tion is quite high not only among MS galaxies (10-11/25, 
~ 40 ± 10%) but also among sub-MS galaxies, which are ei¬ 
ther quenching or already quenched (7/31 ~ 22 ± 7%). One 
would have expected that fading of star formation should be 
accompanied by fading of nuclear activity as well. Why this 
is not the case? One possibility is the following. Galaxies 
may be quenched because at some point they are evacuated 
of virtually all their gas content, or gas supply from the cir- 
cumgalactic medium is somehow suppressed. Yet, this does 
not mean they will remain gas free forever: even in absence 
of gas infall from the environment. The mere gas return 
from dying stars is sufficient to replenish the gas reservoir of 
galaxies over relatively short time scales, with gas accumu¬ 
lating in the potential well. Even if supernovae of Type la 
may provide enough energy to drive galactic winds and out¬ 
flows in less massive galaxies, in the most massive quenched 
galaxies the gas accumulation leads inevitably to a central 
cooling catastroph e and the establis hment of an inflow to¬ 
wards the center dCiotti et al.l Il99ll ) . What happens as a 
result of such a runaway cool ing has been explored w i th hy - 
drodynamical simulati ons by Ciotti & Ostrikcr 


and more recently by iNovak, Ostriker fc Ciottil ( 20121 . and 
references therein). The cooling gas is found to feed the cen¬ 
tral SMBH and its radiative and mechanical feedback is able 
to temporarily revert the inflow into an outflow, possibly 
even restoring a global galactic wind. As the galaxy is evac¬ 
uated of its gas, the AGN subsides and gas starts to accu¬ 
mulate again until the next cooling catastrophe takes place 
and another cycle begins. The duty cycle of this intermittent 
behavior is shorter the higher the rate of mass return from 
dying stars, which decreases with increasing age of a stel¬ 
lar population. Thus, the galaxies with the youngest stellar 
populations among quenched galaxies, have the highest rate 


of mass return, and therefore the shortest AGN duty cycle. 
At least qualitatively, this scenario accounts for the high 
frequency (~ 22%) of AGN among the quenched/quenching 
galaxies in our sample. Actually, their low level of star for¬ 
mation, rather than to quenching, could also be due to mi¬ 
nor star-formation activity that may accompany the central 
cooling catastrophe. This intermitting AGN activity is ac¬ 
tually necessary to keep quenched a galaxy once it has been 
quenched in the first place. 


We point out that a large fraction of X-ray detected 
AGNs in quenched galaxies at 1.5 ^ z ^ 2.5 has been 
also recently uncovered by other authors dTanaka et alJ 


20131: Olsen et al. 

120131:1 Whitaker et al. 20131). In particular 

Olsen et al.l (2015 

) found that the fraction of luminous AGN 


(LA'o. 5 —8fcev > 3 x 10 42 ergs -1 ) among quenched and star¬ 
forming galaxies with M* > 5 x 10 10 Mg are comparable 
(i.e., 19± 9%, and 23±5%, respectively) in the 464.5 arcmin 2 
CDF-S survey, which currently provides the deepest 0.5-8 
keV data. They also found an X-ray excess by stacking the 
X-ray undetected quenched galaxies, and claim that 70%- 
100% of them should host a high- or low-luminosity AGN 
(LA'o. 5 —sfcev ~ 10 42 ergs -1 ). The larger fraction of lumi¬ 
nous X-ray AGNs that we found among star-forming galax¬ 
ies, with respect to that among quenched galaxies, mostly 
depends on the different selection criteria that we used to 
separate the two classes, not only based on UV J rest-frame 
colors, but also on the sSFR (cf. Section [4j. By adopting a 
pure UV J color-criterion one finds that the luminous AGN 
fraction among quenched and star-forming gal axies, are ~ 
26±8%, and 28±9%, respectively, nearly as in I Olsen et al.l 
ll 20131) . 


In our massive sample at 1.4 ^2^2, AGNs 
quenched galaxies only account for ~ 5 ± 5% of the to¬ 
tal X-ray luminosity density. This finding a greees well with 
recent resu lts from iRodighiero et al. l (120151 ) at similar red- 
shift, and iGeorgakakis et al.l ( 2014) at 2 < 1.2, show¬ 
ing that the bulk of the accretion density of the Uni¬ 
verse seems to be associated to normal star-forming (i.e., 
MS) galaxies. For MS galaxies we estimated an average 
log(Lx/SFR) =41.62 ±0.08, which is from 0.6 to 0.8 dex 
higher than the average valu e s foun d by iMullanev et al.l 
d2012al ). and IRodighiero et al.l d2015l ) for similarly massive 
MS galaxies at z ~ 2. This discrepancy comes from the av¬ 
erage X-ray luminosity being derived by combining both de¬ 
tected and stacked X-ray signals, while here we only consider 
X-ray detected objects. A better ag reement is found w i th the 
average L^/SFR ratio derived by IRodighiero et al.l d2015l ) 
for X-ray detected objects in their two higher mass bins 
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(i.e., log(Lx/SFR) = 41.35 ± 0.05, 41.38 ± 0.05, reported 
to a Salpeter IMF, G. Rodighiero private communication). 
The residual discrepancy (~ 0.3 dex) should be attributed 
to the different depths of the X-ray surveys in the GOODS 
and COSMOS fields. 

7 A LESSER IMPACT OF ENVIRONMENT ON 

QUENCHING AT Z = 1.4 — 2 

7.1 The z ~ 1.61 overdensity 

According to the IPeng et al.l (l2010bl 'l phenomenological 
model, environment quenching is not expected to play a 
major role at the redshifts and stellar masses explored 
in this paper. The existence of a diffuse, and well doc¬ 
umented sheet-like over density at z ~ 1.61 on the en- 
tire GOODS-South fiel d (ICastellano et al.ll2007l : iKurk et al.l 

2009; Gia valisco et al.l l201lh offers the opportunity to test 
this expectation by comparing the abundance of quenched 
galaxies, bulge-dominated galaxies, and X-ray detected 
AGNs, within a n d out side the overdensity. We note that 
ICastellano et all (l2007l l isolated a symmetric structure cor¬ 
responding to the highest density peak, centered approxi¬ 
mately at R.A.= 03 ,l :32 m :29.28 s , DEC.= -27°:42':35.99", 
and with an extension of ~ 3 x 3 Mpc comoving, embedded 
in such a wide wall-like or filamentary overdensity. They 
considered as possible members of such cluster-like struc¬ 
ture only the objects within a square of side 2 Ra (where 
Ra = 2.14 Mpc is the Abell radius). 

Here we consider as part of the overdensity all the 
galaxies with spectroscopic redshift in the range 1.6 < 
Zspec ^ 1.62 (19 sources), or photometric redshift in the 
range 1.5 < z p hot < 1-7 (8 sources). Then, we also tested 
the results by only considering objects in the highest den¬ 
sity peak as in ICastellano et alj (l2007h (16 in total, 11 of 
which spectroscopically confirmed). The results are summa¬ 
rized in the following, where the fractions relative to the 
highest density peak are shown in parenthesis. 

• The fraction of quenched galaxies (i.e. galaxies with 
log(sSFR/sSFRMs) < —1) within the overdensity is 59±9% 
(68 ± 12%), which is comparable to the fraction at all other 
redshifts, outside the overdensity, which is 52±9% (50±8%). 

• The fractions of bulge-dominated objects ( B/T > 0.5) 
within, and outside the overdensity are also comparable 
within the errors, being 59 ± 9% (62 ± 12%), and 41 ± 9% 
(45 ± 8%), respectively. They are also comparable to the 
fraction of bulge-dominated galaxies found by Bruce et al. 
(2012) in the same mass and redshift range, in the UDS field 
(55 ± 9% at 1.4 < z < 2, or 65 ± 9% at 1.5 < z < 1.7). 

• Out of our 18 X-ray detected AGN hosts, 12 (6) are in 
the overdensity, and 6 (12) in the field. The relative fractions 
of X-ray detected AGN are ~ 44 ± 9% (37 ± 12%) within 
the overdensity and ~ 20 ± 7%(30 ± 8%) in the field. 

In the overdensity there appears to be a slightly higher 
fraction (at the ~ la level) of quenched, bulge-dominated 
and X-ray detected galaxies, but given that these proper¬ 
ties correlate with each other we cannot claim to have un¬ 
ambiguously detected an environmental effect. Overall, this 
indicates a lesser role for environment quenching with re¬ 
spect to the mass quenching within the studied redshift and 



z 


Figure 10. Redshift distribution of the sample. The diagonal- 
shaded red histogram shows the distribution for spectroscopic 
redshift only, whereas the vertical-shaded distribution refers to 
the addition of galaxies with only photometric redshift 

mass ranges, where no morphology-density relation is yet 
observed. Concerning the X-ray detected AGN, we found 
different results when considering the more diffuse sheet¬ 
like structure, and the highest density peak, since 6 AGNs 
are at a distance of ~ 8' — 10 / (~ 4-5 Mpc, southward ) from 
the cluster center identified by lCastellano et al.l (|2007h . This 
discrepancy is probably due to statistical uncertainties due 
to the small number of objects in our sample. So, our data 
do not allow us to establish whether or not the overden¬ 
sity has an effect on the freque ncy of AGN activity, of the 
kind proposed by some authors llRuderman fe Ebelinsll2()05l : 

iFassbender. Suhada fe Nastasill2012l 'l. 


8 DISCUSSION AND CONCLUSIONS 

This paper is meant to be a step towards shedding new light 
on the physical process(es) responsible for the suppression 
of the star-formation in the most massive galaxies at high- 
redshift, leading to the build up of the high-mass end of 
the quenched population observed in the local universe (i.e., 
the so-called mass quenching). To this purpose, we took ad¬ 
vantage of the high quality WFC3 /HST near-IR imaging 
(CANDELS), as well as of the multi-wavelength data (from 
X-ray to radio) available in the GOODS-S field, to study 
a mass-complete sample of galaxies at 1.4 ^ z ^ 2 with 
M* ^ 10 11 Mq. At this cosmic epoch the SFR density and 
the AGN activity reached their peak and, at the same time, 
the quenching of star-formation started to be very efficient 
in galaxies with mass > 10 11 Mq. Hence, we chose this mass 
and redshift ranges so to deal only with galaxies which are ei¬ 
ther already quenched, or which have initiated the quenching 
process or soon will do so. The purpose is indeed to study the 
interconnection between mass quenching , and AGN, star- 
formation activity, and morphology, especially with the aim 
at identifying transition objects that have left the main se¬ 
quence since a short time. Such transition objects are cer- 
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tainly rare compared to fully star-forming, main sequence 
galaxies or fully quenched ones, hence simple automatic cri¬ 
teria for identifying them are likely to be contaminated by 
these more common neighbors. 

Although many studies have been performed in similar 
mass and redshift ranges, most of them were based on very 
large samples, which while providing statistical significance 
could lead to miss relevant information if automated pro¬ 
cedure are used for the analysis. Instead in this work we 
study galaxies one by one, so to reliably derive their spe¬ 
cific SFR by addressing blending problems in mid-/far-IR 
data, and comparing results from different diagnostics. The 
studied sample consists of 56 galaxies, culled from the K- 
selected D07 optical-to-IRAC catalog, so to include all the 
galaxies with M* ^ 10 11 M© (Salpeter IMF) and spectro¬ 
scopic or photometric redshifts at 1.4 ^ z ^ 2. The opti- 
cal/IRAC data were complemented with Spitzer/ MIPS and 
Rersc/iel/PACS/SPIRE data. The main results can be sum¬ 
marized as follows: 

• By comparing MIPS and high-resolution HST images, 
we identified sources which are blended in Spitzer data, and 
used F160W/WFC3/ HST prior positions to estimate the 
amount of flux associated with each blended galaxy. We 
found that the automated MIPS counterpart association 
based on IRAC priors had failed for 7/56 (~ 12%) sources, 
which had been erroneously considered as MIPS-detected, 
or, just in one case, MIPS-undetected. Based on these re¬ 
sults, we derived the relative fractions of objects which are 
PACS-detected (PACS-d, 18/56), MIPS-detected/PACS- 
undetected (MIPS-d, 10/56), and MIPS/PACS-undetected 
(MIPS-u, 28/56). While ~ 12% may appear a small con¬ 
tamination, it is worth emphasizing that these reclassified 
galaxies represent almost the totality of potential transi¬ 
tion/quenching objects in the sample. 

• We examined the position of our galaxies on the 
M t — SFR plane, relative to the main sequence (MS) of star¬ 
forming galaxies at ( z) = 1.7 (the mean redshift of the sam¬ 
ple), formally including all the galaxies with a sSFR within 
a factor of 4 above and below the average relation (the “MS 
band”). To accurately establish the “degree of quenching” 
of each object, we first derived the SFR as the sum of the 
UV and IR contributions (SFR(IR+UV)), and then cross¬ 
checked the results using optical-to-IRAC SEDs, and two- 
colors BzK and UVJ diagrams. Consistently with the idea 
that the mass quenching process is more efficient in the most 
massive galaxies, we observe a turn down of the average 
SFR at such high-masses. In fact, 31/56 galaxies (~ 55%) 
are located below the lower boundary of the “MS band”, 
and 2/3 of the MS galaxies lie below the average M*—SFR 
relation. All the galaxies below the “MS band” turned out 
to have a specific SFR more than 10 x lower with respect to 
the average relation (log(sSFR/sSFRMs) < —1), as well as 
properties typical of already quenched objects, such as SED, 
colors and spectra, when available). This suggests that the 
mass quenching may be a very fast process in massive galax¬ 
ies, being able to shut down the star-formation in relatively 
short timescales (<< 1 Gyr), though a much bigger sam¬ 
ple should be analyzed in a similar fashion to achieve firmer 
conclusions in this respect. 

• We fitted the surface brightness of galaxies with both 
a single-Sersic function and with a double Sersic, in the lat¬ 


ter case to model bulge+disk profiles. We found that the 
bulge to total light ratio ( B/T ), used as a proxy for the 
importance of the bulge component in a galaxy, provides 
results in general agreement with those obtained from the 
single-Sersic index (n). In the full sample ~68% of the ob¬ 
jects have a relevant bulge component (B/T > 0.3), and ~ 
50% are bulge-dominated ( B/T > 0.5). These findings are in 
agreement with recent results from Bruce et al. (2012, 2014), 
showing that while at z > 2 galaxies with M* > 10 11 M@ 
are mostly disk-dominated, at 1 < z < 2 they are primar¬ 
ily bulge+disk systems. In particular, ~84% (26/31) of the 
quenched galaxies have a relevant bulge component, and 
61% (19/31) are bulge-dominated. On the other hand, about 
half of the MS galaxies have B/T > 0.3, and ~36% (9/25) 
are bulge-dominated. The average (B/T) seems to increase 
within the MS, with decreasing sSFR, and then remain al¬ 
most constant below it. Although we do not have enough 
statistics to establish the solidity of this finding, this fits 
well with recent observational results showing that the bulge 
grow th is directly related to the quenching of star-formation 
(e.g. iGenzel et all l2014al : lLang et al.l 120141 ; iTacchella et al.l 
l2014h . 

• We identified 17 secure and one candidate AGN hosts 
in our sample (i.e.~ 30%), both among MS (11/25, ~ 44%) 
and among quenched (7/31, ~ 22%) galaxies, by comparing 
the SFR. derived from Chandra/X- ray and VLA/radio data 
with those inferred from IR+UV, and SED fitting. Although 
the fraction of AGNs among quenched galaxies is quite high, 
their contribution to the X-ray luminosity density is small, 
representing only 5 ± 5% in our sample of M* ^ 10 11 M@ 
at z = 1.4 — 2 , in a greement with recent results from 
iRodighiero et all d2015l l. In general, we found that the pres¬ 
ence of the AGN correlates with a steep central profile, with 
Sersic n ^ 2 — 2.5, and B/T > 0.3 —0.5, not only in quenched 
galaxies, but also in almost all the MS galaxies. This is con¬ 
sistent with the idea that the galaxy bulge and the central 
SMBH grow together since early epochs, as expected based 
on the bulge mass-SMBH relation observed in the local Uni¬ 
verse. Moreover, it seems that AGNs are much more fre¬ 
quent in galaxies with —0.6 <log(sSFR/sSFRMs) < 0 than 
in those with higher specific SFR (0 <log(sSFR/sSFRMs) < 
0.6), suggesting that the AGN is somehow related to the de¬ 
cline of the SFR. This findings provide some support to the 
notion that AGN feedback may be crucial in shutting down 
the star-formation in the most massive galaxies, as assumed 
in some theoretical and semi-analytical models. 

• With respect to the local ETGs, the majority of the 
quenched galaxies (28/31, ~90%) are undersized (by a fac¬ 
tor ~ 2.5, on average), while the majority of the MS galaxies 
(20/25, ~ 80%) are comparable to them in size, in agreement 
with what found in previous works. All the 5/25 compact 
MS galaxies are similar to the “compact star-forming galax¬ 
ies” (CSFG) studied by Barro et al. (2013, 2014), showing 
high central concentration (n > 2.5, and B/T > 0.5) and 
a sSFR slightly lower with respect to the average MS rela¬ 
tion. Moreover, 4 out of 5 are AGN hosts, although in two 
objects the AGN dominates the photometric light, possibly 
affecting the derived galaxy size. All these characteristics 
suggest that these objects will be quenched soon, originat¬ 
ing compact quenched remnants, as suggested by Barro et 
al. 

• In our sample 27 out of 56 galaxies are included in the 













18 C. Mancini et al. 


well-known overdensity at z = 1.61 in the GOODS-S field, 
hence we compared the relative abundances of quenched 
galaxies, bulge-dominated ( B/T 0.5) galaxies, and X-ray 
AGNs, within and outside the overdensity. This test showed 
that no appreciable difference exists in the relative fractions 
of these three sub-populations in the different environments. 
The result does not change significantly when the test is re¬ 
peated by considering as part of the overdensity only the 
objects included in the highest-density region (~ 3 x 3 Mpc) 
identified by Castellano et al. (2007). This suggests that the 
environment should have a minor (if not negligible) role in 
triggering the quenching of the star formation in the most 
massive galaxies at 2 = 1.4 — 2, as predicted by the Peng et 
al. (2010, 2012) phenomenological model. 

Perhaps the most relevant lesson learned from this pilot 
investigation is that careful object-by-object analysis is nec¬ 
essary to properly identify candidate galaxies on their way 
to be quenched, whose number and structure can be im¬ 
portant for the identification of the physical mechanism(s) 
responsible for the quenching of star formation in massive 
high-redshift galaxies. 
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Table 1. Massive Sample: multi-wavelength and spectroscopic information from GOODS/CDFS/CANDELS Surveys. 


RA 

[deg] 

DEC 

[deg] 

z 

zqtD 

Reft 2 ) 

^24 

[hJy] 

-^24, new 

[/jJy] 

Tl.4GHz (3) 

My] 

^Hard —X ( 4 ) 
[10-16 _«|_] 

L cm-* 1 s J 

ix (4) 

[10 42 2£S] 

class(X)/(o)( 5 ) 

sub-sample 

53.061039 

-27.69352 

1.44 

0 

b, CaHK 

< 17.22 






MIPS-u 

53.173473 

-27.697678 

1.427 

1 

b, CaHK-Mgll-Mgl 

155.19 ± 10.64 

< 17.66 





MIPS-d—m 

53.1736946 

-27.6981106 

1.59 

0 


< 31.92 

< 21.0l( 8 ) 





MIPS-u 

53.1411934 

-27.7011375 

1.604 

1 

a 

66.1 i 6.02 



< 4.2 

11.0i3.2 

AGN 

MIPS-d 

53.1763992 

-27.7011547 

1.981 

1 

a 

24.89 i 6.73 






MIPS-d 

53.104889 

-27.705238 

1.617 

1 

d 

564.94 i 19.92 

423.76 ± 11.71 

88.5 ± 14.6 

27.Oil.7 

77.0i2.1 

QSO-l/BLAGN 

PACS-d 

53.1046181 

-27.7054272 

1.609 

1 

b 

< 59.76 






MIPS-u 

53.0620193 

-27.706501 

1.89 

0 


< 28.23 

204.38 ± 9.41 


3.8il.6 

8.9i3.3 

AGN 

MIPS-u—>d 

53.1085548 

-27.7101555 

1.609 

1 

a 

< 18.75( 6 > 






MIPS-u 

53.1830406 

-27.7090015 

1.76 

0 


23.02 i 7.37 






MIPS-d 

53.1496506 

-27.7113838 

1.612 

1 

a 

30.7 ± 6.0 

< 41.55 





MIPS-d—m 

53.179775 

-27.711674 

1.416 

1 

a 

< 19.34 






MIPS-u 

53.089447 

-27.7115822 

1.72 

0 


242.9 ± 14.82 






PACS-d 

53.1164207 

-27.7127037 

1.61 

1 

a 

< 19.78( 6 > 






MIPS-u 
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-27.713728 

1.609 

1 

a 

'lO 

00 

V 
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5.7il.5 

AGN 

MIPS-u 
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-27.7162361 

1.614 

1 

a 

< 12.33( 6 ) 






MIPS-u 
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1.66 

0 
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< 19.18(8) 





MIPS-d—m 

53.130501 
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1.415 

1 

a 

< 25.76( 6 ) 






MIPS-u 

53.188858 
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0 
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MIPS-d 
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1 

c 

289.09 ± 16.27 
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240.0i4.0 

QSO-l/BLAGN 

PACS-d 
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1 

a 
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0 
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2 
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2 

c 
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MIPS-d—m 
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1 

c 
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AGN 

MIPS-d—m 
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1 

c 

254.29 ± 15.24 






PACS-d 

53.0521927 

-27.7747669 

1.605 

1 

c 

< 19.06( 7 ) 



2.0i0.8 

4.9il.3 

AGN 

MIPS-u 

53.125061 

-27.7907791 
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1 
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MIPS-d 

53.0446396 

-27.8054123 

1.61 

0 
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MIPS-u 

53.101093 

-27.808559 

1.97 

0 


116.99 A 8.68 






PACS-d 
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1.413 

1 

b, [OII]-D4000 
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<1.0 


PACS-d 
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0 
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MIPS-d 

53.1572495 

-27.8335114 

1.619 

1 
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MIPS-d—m 
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1 
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PACS-d 
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-27.834509 
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1 

b, [OII]-MgII 

< 31.10 






MIPS-u 

53.1271667 

-27.8345394 

1.81 

0 


< 15.28(8) 






MIPS-u 

53.1314774 

-27.8412781 

1.613 

1 

b, OII-Mgll 

191.47 ± 12.52 
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1.2i0.5 

2.1i0.7 

AGN 
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53.168438 

-27.8618145 

1.77 

0 


< 32.90(®) 






MIPS-u 

6898 

53.1573715 

-27.8700867 

1.603 

1 

e 

1085.15 ± 25.54 


112.7 ± 12.7 

98.0±2.8 

250.0±3.7 

QSO-l/HEX 

PACS-d 

7077 

53.2080154 

-27.8743706 

1.609 

1 

b 

104.89 ± 8.13 






MIPS-d 

7116 

53.1824493 

-27.8753071 

1.48 

0 


163.92 ± 11.1 






PACS-d 

7491 

53.148407 

-27.885511 

1.604 

0 


< 12.27< 6 ) 






MIPS-u 

7617 

53.099564 

-27.890564 

1.57 

0 


< 18.2l(®) 






MIPS-u 

7923 

53.079632 

-27.901369 

1.4 

0 


< 50.40 

< 17.12< 6 > 





MIPS-u 

8121 

53.2112389 

-27.9082909 

1.81 

0 


< 15.7l(®) 



< 4.8 

4.8±1.6 

AGN 

MIPS-u 

8569 

53.121048 

-27.9285126 

1.469 

1 


182.92 ± 12.09 






MIPS-d 


(^Redshift quality flag: O=photometric, l=spectroscopic, secure, 2=spectroscopic, likely (in agreement within 68% confidence limits with the photometric redshift). 

References for spectroscopic redshifts: (a)GMASS (KurK et al. 2012), (b)FORS2 (v3.0, Vanzella et al. 2008), (c)K20 (Mignoli et al. 2005), (d)X-ray CDFS (Szokoly 
et al. 2004), (e)VIMOS (Popesso et al. 2009; Balestra et al. 2009). 

( 3 )VLA Radio integrated flux Data from 1 Keller mann et alj (120081) , ATCA survey, [Norris et al.l (120061 ), and I Afonso et all d2006l) 

( 4 ^Data from 4 Ms Chandra X-ray catalog dXue et al. I fioTil) : Flux Hard-X band(2-8 kev), and total integrated Full-band (0.5-8 kev) luminosity. 

Objects identified as AGN, based on X-ray/Radio data as explained in Section l3.3l X-ray/Optical classification from Szokoly et al. (2004) is also reported, 
when available: class(X)=X-ray class (AGN-1,2, or QSO-1,2); class(o)=classification only based on optical spectra (HEX= object with unresolved emission lines and 
exhibiting high ionization lines, or emission line ratio indicating AGN activity. These objects are dominantly optical type-2 AGNs or QSOs; BLAGN= objects with 
emission lines broader than 2000 km/s. This classification implies an optical type-1 AGN or QSO). 

( 6 )2ct upper limits in MIPS for objects with a S/N < 1 (3cr upper limits are used for the remaining objects). 

(^Marginal MIPS detection, i.e., 2.5<S/N<3, treated as 3 a upper limit. 

( 8 ) This galaxy has been classified as MIPS-u, being very likely a lens system, whose MIPS flux (Flux(24 fim) n ew =152.12 ± 7.0), derives from a background source 
likely at z ~ 2 — 2.5 (see Appendix 1X1). 
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Table 2. SFR, Stellar Mass, and SED fitting parameters (BC03, models, Salpeter IMF). 


ID 

SFR(IR+UV) 

Af* (D04) 

AM SED) 

Age 

A V 

SFR(SED) 

T 

T min T max 


[M @ /y r ] 

[-Mg] 

[M 0 ] 

[Gyr] 

[mag] 

[M@/y r] 

[Gyr] 

[Gyr] 

282 

<12.89 

11.18 

n 17 +o.oi 

1 ' —0.02 

2 10+ 000 
ZiU -0.10 

0 10+ 0 - 05 
ulu -0.05 

0 62 +°- 22 

0-0^ -o.03 

0.3 

0.3-0.3 

427 

<12.95 

11.26 

11 11+ 002 
11,1 —0.01 


0 50"*”°' 05 

U.OU_o ig 

9, +0.07 
-0.02 

0.1 

0.1-0.1 

428 

<26.31 

11.04 

10 99"*" 0 ' 03 
xu.yy _ Q 01 

3.7515“ 

1 io+ 0 - 15 

1U -0.15 

33.1115° 3 5 2 

30.0 

3.0-1000 

552 

56.23+5.09 

11.17 

u-«±S;i 8 i 

0 72+ 1 - 18 
u - ‘ z - 0.08 

1 20 - * -0 ' 15 
iZU _0.20 

59 qq +20.61 
uy.yy —12.36 

0.3 

0.3-1 

557 

18.08+4.89 

11.01 

n-i8±g;“? 

1 61+ 007 
1DJ --o.n 

0 ^O - * -0 ' 05 
U.OU_o 05 

3 19 + 013 
—0.62 

0.3 

0.3-0.3 

686 

332.14+9.39 

11.04 

11 04"*" 0 ' 01 
—o 00 

0 45+ 0 00 
u -^°_0.01 

0 OO - * -0 ' 03 
u.uu-o 00 

0 00 +0 00 

u.uu -o.oo 

0.0 

0.0-0 

687 

< 47.23 

11.30 

u-34±S;S? 

2 io+ 010 
2iU -0.08 

0 30 - * -0 ' 05 

U - 3U _0.05 

0 90 + 0,0 ^ 

u.yu _o 24 

0.3 

0.3-0.3 

720 

132.89+6.46 

11.06 

n-osto.oi 

0 64"* -0 ' 57 
-0.00 

0 60"*” 0 01 
u ' ou —0.30 

9 co +0.36 

Z 0,5 -0.03 

0.1 

0.1-0.2 

848 

<14.65 

11.18 

n n+o.04 
—0.03 

1 oq+0.16 

-0.00 

0 40"*"°'° 5 

U - 4U _0.05 

0 00 +0-00 
u.uu -o.oo 

0.0 

0.0-0.1 

856 

15.09+4.90 

11.57 

11 52+ 0 01 

n.oz_o.oo 

i oa +0.00 

i - 28 -o.oi 

0 oo+ 002 
u.uu_o 00 

0 00 + 0 - 00 

u.uu -o.oo 

0.0 

0.0-0 

880 

<32.88 

11.23 

11 21"*" 0 ' 02 

11 zl _o.01 

1 14+0.00 
— 0.00 

1 65"*" 0 ’ 05 
i.DO-o 04 

16 06 +°- 64 

10.uo -0.50 

0.3 

0.3-0.3 

887 

<13.94 

11.23 

11 05~*" 0 ' 02 
ll.UO_o oi 

1 9Q + O.OO 

J-^O-O.oo 

0 io+ 0 05 
ulu -0.02 

0 0 
0 0 
d d 
+ 1 

O 

q 

6 

0.0 

0.0-0 

903 

164.43=1=36.70 

11.08 

11 f)6~*~ 0 ' 05 

ii.uo_q 03 

n i a+o.ii 
Ui -0.02 

9 ye+0.05 
/o -0.05 

269.90 1 IIH 

0.1 

0.1-0.3 

947 

<15.44 

11.06 

11 10+0- 01 

1 90"*"°' 10 
i.yu_o 10 

0 30 - * -0 ' 05 

U.8U-0 10 

1-00 15;55 

0.3 

0.3-0.3 

986 

<14.65 

11.42 

11 33+ 001 
n.oo_o 02 

1 9O+0.00 

J-^O-O.oo 

n qc + 0.03 
U.30_o 05 

0 00 +0-00 

U.UU _Q QQ 

0.0 

0.0-0 

1084 

<9.58 

11.21 

11 22~*~ 0 ' 02 
n.zz_o oi 

2 10+ 0 - 10 
ZiU -0.10 

0 15+ 0 - 10 
ulo -0.05 

0.69 15“ 

0.3 

0.3-0.3 

1187 

<14.38 

11.24 

11 26 - * -0 ' 0 '*" 

11. zo — o.Ol 

1 61+ 0 - 04 
101 — 0.00 

0 45+0- 02 
u -^°_0.06 

3.80 15'“ 

0.3 

0.3-0.3 

1272 

<18.00 

11.06 

10 9i+0- 01 
iu.yi _q 02 

1 9Q + O.OO 

J-^O-O.oo 

0 40"* -0 ' 04 

u -^ u _0.05 

0 00 +0-00 

U.UU _Q QQ 

0.0 

0.0-0 

1510 

238.82=1=13.20 

11.08 

11 oi"*' 002 

0 23"*” 0 ' 03 
U - Zl3 _0.02 

i-yu_o .10 

141.22 155“ 

0.1 

0.1-0.1 

1906 

123.32=1=24.14 

11.04 

10 74"*" 0-01 
iu. '^_ 0 .00 

0 io"*" 0-00 
ulu -0.00 

1 65"*” 0 ' 02 
i.oo_o 00 

62i - s3 mil 

1000.0 

1.0-1000 

2450 

<10.44 

11.02 

11 14+ 001 

1 —0.02 

i-6<°:oo 

0 55"*" 0 ' 02 
u.oo-o 10 

9 07 +0.06 

2 - 87 -0.75 

0.3 

0.3-0.3 

2470 

133.70=1=43.62 

11.02 

11 31+004 

ll.cJl_o go 

3.5015;“ 

9 9C+O.8O 
-0.20 

41 31 + 8114 

41.8± _ 14 64 

3.0 

0.1-1000 

2940 

<12.10 

11.22 

Ii.30±g;g* 

1 tjo+0.03 

128 -0.14 

0 60"*”°' 25 
u.ou_o 04 

12 67 + 9 - 25 

LZ.Ot _Q 43 

0.3 

0.3-0.3 

3066 

93.05+27.05 

11.12 

i0.84+°;°2 

0 io"*" 0-01 
ulu -0.00 

3 00+ 0 - 00 
3.uu_o 05 

770 96 +1^-06 
( /'u.yo 10 5.34 

1000.0 

0.3-1000 

3230 

<32.67 

11.09 

H-°°±0.02 

-1 9Q + O.OO 

J-^O-O.oo 

0 15+0-05 
U-lo -0.05 

0.001555 

0.0 

0.0-0 

3231 

<37.14 

11.65 

ii-68±°;°? 

1 on"* -0-06 
i.yu_0 00 

0 55+0-05 

U.OO-o 06 

3.801515 

0.3 

0.3-0.3 

3236 

114.20=1=16.96 

11.03 

10 98”*~ 0 ' 01 

iu.y<5_ 0 Q4 

0 20~*~°' 01 

U '^-0.02 

1 70+°- 10 
-*-• ‘ u —0.02 

176 47 + 45 - 19 

i/o.4f 10 98 

0.1 

0.1-0.1 

3258 

<14.91 

11.31 

11 qo+0.01 

11.38_o 02 

2 00"*" 0-10 
^.uu_o 10 

0 35"*" 0 ' 05 

U.80_o IQ 

. 00 +0.49 
-0.41 

0.3 

0.0-0.3 

3853 

23.44=1=6.11 

11.31 

H-SO+o'qj 

1 14+O.OO 
-0.04 

1 05+ 0 - 07 

l.UO-o 05 

20 14 + 4 05 

zu.14 _ 1 QQ 

0.3 

0.3-0.3 

4310 

<10.49 

11.03 

10-95^:o2 

■1 90+O.OO 

i - 28 -0.00 

0 40"*'®' 33 

-0.05 

0 00 +0 00 
u,uu -0.00 

0.0 

0.0-0 

4406 

137.21+63.97 

11.21 

11 24~*" 0-02 

1 14+ 1-31 
— 0.00 

0 55+O.IO 
Z.OO-o 05 

1 7 38 +0.77 
!'• 38 -0.39 

0.3 

0.3-0.3 

4705 

199.93=1=5.83 

11.01 

10 91+ 0 03 

iu.yi_o 03 

n 90+O.O3 

9 or\+0.10 

2.8U_o 05 

119 53 +30.33 

11Z.88 23.24 

0.1 

0.1-0.1 

5149 

<12.22 

11.10 

11 14+OOO 
—0.02 

1 on+0-07 

1.80-0 12 

0 40"*" 0 ' 10 

u -^ u _0.05 

I BS toll 

0.3 

0.3-0.3 

5415 

39.38+5.00 

11.05 

11 12+0-00 
lllz -0.03 

o-8+5:S 5 9 

1 oo+ 0 - 10 

i.uu-o os 

40.98 1“ 7 5 6 

0.3 

0.3-0.3 

5503 

426.87=1=9.27 

11.06 

1 0 9^+0 05 
lU.yD_ 0 02 

0.1815;“ 

1 80+®'® 3 
l-od_0.15 

211 19 + 2 334 
Zll.iy 72.04 

0.1 

0.1-0.1 

5509 

<13.54 

11.32 

11 54+ 0 03 

0 7C+O.OO 

° —0.25 

. HC -+0.10 
1.8&_o 05 

11 +2-99 

1146 -0.34 

1.0 

1.0-1 

5530 

<27.52 

11.10 

10 94"*" 0-01 
iu.y4_o Q1 

1 70+°- 10 
' -0.02 

0 30"*”°' 05 
O-5O-0.10 

, oe +0.13 
-0.38 

0.3 

0.3-0.3 

5534 

375.88=1=10.78 

11.23 

11 23~*~ 0 ' 02 
ll.Z3_o oi 

n 09+O.OO 
u - l3Z _0.03 

2 45+ 0 - 15 
^• 4O _0.03 

88.29 15 9 13 3 

0.1 

0.1-0.1 

5556 

<24.83 

11.08 

11 03"*" 0 ' 02 

11. Uo — o.Ol 

n si +o.oo 
u -°-o.oo 

0 45 - * -0 ' 04 
o-^°_o.03 

0 43 + 0 - 02 

-o.oi 

0.1 

0.1-0.1 

5564 

<9.34 

11.22 

11 34+ 0 01 
11.84-0 02 

2 00"*”°' 10 
Z.UU-o 07 

0 30" 1 " 0 ' 07 

o.8U_ 0 1Q 

-1 97 +0.38 
-0.38 

0.3 

0.3-0.3 

5860 

141.95=1=15.23 

11.40 

11 qq+0.02 

11.80 —ooi 

0 29"*" 0 ' 00 

u z -0.03 

2 45+ 0 - 15 
z - 4O _0.06 

158.6015155 

0.1 

0.1-0.1 

5974 

511.86=1=24.04 

11.01 

11 20~*~ 0 ' 03 
ii.zU — o 09 

0 64~*~ 008 

U- ° -0.35 

2 20 - * -0 ' 25 
zzu - 0.10 

90 06 + 22 - 89 

yu.uo _2i.o8 

0.3 

0.1-0.3 

6071 

185.39=1=37.17 

11.29 

11 10+0-05 

ii.iU-Q 06 

0 20”*~ 0 ' 12 
u -^ u _0.07 

9 qc + 0.05 
2.8t)_ 0 05 

51 q 4« +141.29 
913.40 _ 138 .94 

0.3 

0.1-1000 

6193 

126.22=1=8.65 

11.32 

ii-i9i°o: 0 o°o 

0 64~*~ 000 
u o -0.01 

0 60"*” 0 ' 00 
U.OU-o 01 

0 no +0.00 

- 3 - 420 -0.02 

0.1 

0.1-0.1 

6224 

74.75=1=11.88 

11.02 

10-68tS'.o? 

0 20~*~°’ 08 
u ‘^ u -0.06 

9 oc+0.05 
Z.80_o IQ 

280.64 155;“ 

1000.0 

0.3-1000 

6352 

91.69=1=21.73 

11.03 

10 94+ 0-03 
iu.y4_ 0 oo 

0 32”*” 0 ' 04 

U " 3Z —0.03 

1 15+ 0 - 15 

1 i0 —0.09 

45.71155.55 

0.1 

0.1-0.1 

6572 

265.86=1=77.08 

11.08 

11 50+ 005 

ii.oU — 0 og 

3 00+ 0 0 ° 
3.uu_i 20 

2 00+ 0 15 
z.ou_ 0 20 

132.45 tiUl 

30.0 

3.0-1000 

6647 

<20.92 

11.10 

11 14+0-01 

1 —0.02 

1 oq+0.04 

i 2 -0.00 

0 60"*” 0 ' 00 
u.uu-o 05 

0 00 +0-00 
u.uu _Q 00 

0.0 

0.0-0.0 

6898 

90.28=1=16.12 

11.41 

11 56 +0 05 
11.00 —o 09 

1 43+ 0 - 47 

J--^ —0.79 

9 9.+O.20 
z,zo -0.15 

305.12 l“ 8 6 f 

10.0 

1.0-1000 

7077 

83.23=1=6.59 

11.22 

11 17+ 002 
111 ' -0.00 

0 36”*” 0 ' 00 

U.8D-0 04 

1 65+ 0 - 15 
1.OO-0 01 

51.81 15B 5 “ 

0.1 

0.1-0.1 

7116 

90.94=1=31.19 

11.13 

11 oq+0-01 

ll.zy _o 04 

4 00"* -0 ' 00 

^' uu — 0.75 

1 7C+0.05 
' -0.10 

66 21 + 6 - 91 
00.21 _ 10 Qg 

1000.0 

10.0-1000 

7491 

<9.68 

11.09 

10 98"*" 0 ' 01 
iu.yO —o Q 2 

-1 90+O.OO 

J-^O-O.OO 

0 1 5+0-03 
u - lo -0.05 

0 00 +0 00 

u.uu _Q qq 

0.0 

0.0-0.0 

7617 

< 14.47 

11.03 

10 97+O O 3 
iu.y f _o 03 

-1 9 q-|- 0.23 

J-^O-O.OO 

0 OO"*" 0 00 
U'UU — 0.00 

0 00 +0-00 
u.uu _Q qq 

0.0 

0.0-0.0 

7923 

<11.91 

11.11 

11 01 - * -0 ' 00 
ii.ui-Q 02 

1 90 - * -0 ' 01 
i.yu_o 10 

0 OO"*" 0 05 
U.UU-O 00 

0.83 15“ 

0.3 

0.3-0.3 

8121 

<10.25 

11.37 

11 51+ 002 
ii.oi-Q oi 

1 90"*”°' 10 
i.»u_0 03 

0 25"*" 0 ' 04 
U -^°_0.10 

2 62 +°- 12 
-0.76 

0.3 

0.3-0.3 

8569 

140.69=1=9.35 

11.33 

11 14+0-01 

n 90+O.O3 

o-^^-o.oi 

, 7 c+0.02 

± - /o -0.10 

191 17 + 1401 

-44.35 

0.1 

0.1-0.1 


Note. — The second column shows the total SFR(IR+UV), where the IR contribution is derived from the IR SEDs, 24 /^m/MIPS 
flux, or 24 /i,m/MIPS upper limit, for PACS-d, MIPS-d, and MIPS-u galaxies, respectively, as detailed in Section |3.2| Af*(D04), 
and Af* (SED) are stellar masses derived as in D04, and from the best-fit SED, respectively. Other best-fit SED quantities follow; 
in the order: mean stellar-age, reddening (Ay), SFR, declining t, and r 90% confidence interval. All the uncertainties on the 
SED-derived quantities are within the 90% confidence level, and inferred from the % test. Typical uncertainties on Af* (D04) 
are around 40% (cf., D04). 
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Table 3. Morphological parameters from the GALFIT 2D SB fitting on the WFC3/fl60W (if-band) images. 


ID 

mag 

n 

^e,circ 

B/T 

PSF/T (sS+PSF) 

/ rrT sSFR 

^e,circ 

Luy sSFR(MS ) 

K 0 


[AB] 


[kpc] 




282 

21.31 ± 0.01 

4.68 ± 0.08 

1.63 ± 0.02 

0.83 ± 0.01 

< 0.1 

2 45 + 0 - 16 
^• 4£) -0.02 

0.41 

427 

20.91 ± 0.01 

3.73 ± 0.04 

2.4 ± 0.02 

0.47 ± 0.02 

0.15 ± 0.006 

-2 20 + 0 - 02 
z _o.01 

0.55 

428 

21.65 ± 0.02 

3.8± 0.08 

5.8± 0.17 

0.2±0.01 

0.13 ± 0.01 

-0 6 + 0 08 
u -° -0.06 

1.78 

552 

21.18 ± 0.01 

4.01 ± 0.09 

3.66 ± 0.08 

0.54 ± 0.04 

< 0.1 

-0 47 + 0 - 04 
-0.04 

0.95 

557 

22.28 ± 0.01 

2.26 ± 0.04 

2.07 ± 0.02 

0.65 ± 0.03 

< 0.1 

1 59 + 0 02 

_q.08 

0.66 

686 

21.73 ± 0.02 

4.02 ± 0.19 

2.55 ± 0.05 

0.85 ± 0.04 

0.22 ± 0.01 

0 45 + 0 - 03 
-0.03 

0.78 

687 

21.37 ± 0.01 

2.66 ± 0.02 

1.8 ± 0.0 

0.62 ± 0.02 

< 0.1 

o oo +0.02 
-2.38 —0.12 

0.39 

720 

21.72 ± 0.09 

3.39 ± 0.42 

4.71 ± 0.73 

0.32 ± 0.05 

< 0.1 

-0 01 + 0 - 02 

UU1 -0.02 

1.41 

848 

21.58 ± 0.01 

3.23 ± 0.05 

0.88 ± 0.01 

0.83 ± 0.02 

< 0.1 

<-3.0 

0.23 

856 

20.45 ± 0.01 

5.31 ± 0.04 

2.85 ± 0.02 

0.86 ± 0.02 

0.10 ± 0.01 

<-3.0 

0.44 

880 

21.91 ± 0.01 

4.8 ± 0.16 

1.59 ± 0.03 

1.0 ± 0.0 

< 0.1 

-1 06 + 0 02 
i.uo _ 0 oi 

0.38 

887 

21.04 ± 0.01 

5.9 ± 0.08 

2.38 ± 0.03 

0.76 ± 0.02 

< 0.1 

<-3.0 

0.57 

903 

22.18 ± 0.01 

0.88 ± 0.01 

4.36 ± 0.03 

0.0 ± 0.0 

< 0.1 

0 07 + 0 - 10 
u-u ' -0.10 

1.28 

947 

21.83 ± 0.01 

2.69 ± 0.05 

0.7 ± 0.0 

0.71 ± 0.03 

< 0.1 

9 io +U. 15 
-z.io _ 0 13 

0.21 

986 

20.96 ± 0.01 

6.27 ± 0.09 

1.89 ± 0.02 

0.79 ± 0.02 

< 0.1 

<-3.0 

0.35 

1084 

21.63 ± 0.01 

2.59 ± 0.04 

1.13 ± 0.01 

0.0 ± 0.0 

0.24 ± 0.01 

.2 41 + 019 
-0.12 

0.28 

1187 

21.5 ± 0.01 

2.73 ± 0.02 

2.11 ± 0.01 

0.0 ± 0.0 

0.18 ± 0.01 

-1 94 + 0 01 
-0.14 

0.50 

1272 

21.74 ± 0.01 

2.54 ± 0.03 

2.02 ± 0.02 

0.81 ± 0.02 

< 0.1 

<-3.0 

0.61 

1510 

21.61 ± 0.01 

1.01 ± 0.01 

6.94 ± 0.08 

0.0 ± 0.0 

0.10 ± 0.05 

0 22 + 0 - 02 
-0.02 

2.02 

1906 

21.39 ± 0.01 

5.01 ± 0.05 

1.68 ± 0.01 

1.0 ± 0.0 

< 0.1 

-0 03 + 0 09 
_o.09 

0.52 

2450 

22.12 ± 0.01 

1.9 ± 0.03 
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Note. — The total magnitude, Sersic index, and circularized effective radius are from the single-Sersic fit (see Section |5.U . 
For each object we also report the Bulge to Total flux ratio (B/T) from the composite Bulge+Disk model, discussed in 
Section |5.2| and the fractional contribution to the total flux of a hypothetical central point-source added to the single-Sersic 
fit, PSF/T(sS+PSF), when it is > 10%. The last two quantities are the distance from the MS in log-scale, as shown in 
Figures |5|8| and and the distance from the M* — size relation of local ETGs, shown in Figures^] and fSI 
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APPENDIX A: MIPS DEBLENDING 

In this appendix we detail the procedure used to derive the 
MIPS flux (and SFR) estimates for blended sources in our 
sample, whose results are shown in Section [4] and Figure [2] 
Since the 24 ^im/MIPS FWHM (= 5"9) is substantially 
larger than the size of our galaxies, each blended system 
includes all the objects closer than ~ 3” (HWHM) to the 
centroid of the MIPS emission. We used GALFIT to decom¬ 
pose the blended sources by fitting PSF models on the fixed 
F160W/WFC3 positions, leaving the galaxy magnitude as 
the only free parameter to be recovered. To avoid to have 
residual flux in the neighborhood, the sources at a distance 
3 — 4" < r < 20 " from the galaxy(ies) of interest were also 
fit together, but fixing all the parameters (i.e., position and 
magnitudes) to the known values from the available MIPS 
catalog. In the following we discuss the results derived for 


25 


each object, or pair (when both the galaxies in a blended 
pair are included in our sample). We considered as MIPS- 
detected (to be included in MIPS-d subsample) objects with 
SNR> 3, and treated the others as upper limits (at 2a or 
3a, for SNR< 1 and ^ 1, respectively). 

• Objects #427 and #428. In the MIPS catalog based 
on IRAC positions the MIPS counterpart had been wrongly 
assigned to #427, which is not compatible with such high 
MIPS flux, being spectroscopically classified as Early-Type 
galaxy as proved by the CaHK,MgII,MgI absorption lines 
(GOODS) and lack of the O[II] emission line (expected 
for a star-forming galaxy in the (VLT/FORS2) observed 
wavelength-range dVanzella et alj|2008l h Moreover, the re¬ 
sults of PSF fitting on the image confirm that object #427 
has a low 24 /im flux (F24< 26 /rJy), while most of the 
MIPS flux should come from object #428. However, as de¬ 
tailed in Section ona this is one of the most peculiar case of 
blending in our sample, being consistent with a double lens 
system. In this perspective the MIPS flux should be associ¬ 
ated to a lensed background star-forming galaxy, probably 
at 2 ~ 2 — 3. 

• Objects #686 and #687. These galaxies constitute 
a gravitationally bound system, as proved by their spec¬ 
troscopic redshifts and mutual distance (Appendix [B]). The 
PSF fitting on the J7ST/WFC3 positions suggested that 
most of the MIPS emission comes from the AGN host #686, 
although a small percentage of the flux may be associated 
to the companion (#687). However, since both the red col¬ 
ors of #687, and the SED fitting results, suggest old age 
(~ 2 Gyr), and very low reddening {Ay = 0.3 mag), and 
SFR (~ 1 Mg/yr), we considered this object as MIPS- 
u(cf. Table EJ. We notice that, also if considered as MIPS- 
detection, #687 would remain a “border-line” galaxy more 
than 4x below the MS, since the MIPS flux derived from 
PSF fitting (#24 = 65.2 ± 13.8 /rjy) is very close to the 
quoted 3a upper limit (#24 < 59.76 AtJy). 

• Object #720. This galaxy has a near-IR/optical 
bright and massive (Af* ~ 5 x 10 11 Mg) neighbor (i.e., ob¬ 
ject #721), not included in our sample since its photometric 
redshift is slightly lower than the cut used in our sample- 
selection, i.e., Zphot = 1.36. In the GOODS catalog #721 had 
been identified as the only counterpart of the bright MIPS 
source. On the contrary, our deblending procedure suggested 
that most of the MIPS flux should be associated to object 
#720 (#24 = 204.4 ±9.4 /rJy), with a smaller fraction of the 
flux (#24 = 83.5 ± 10 /rJy) coming from #721. The presence 
of peculiar tidal tail seems to suggest that this could be also 
a merging system, but this hypothesis is not supported by 
spectroscopic information. 

• Object #880. The detection from the GOODS cata¬ 
log (SNR> 3) was replaced with an upper limit, since the 
MIPS source seems to be the result of two blended marginal 
detection, as also suggested by the elongated shape of the 
source in the MIPS image. 

• Object #1187. As for object #880, this MIPS detec¬ 
tion appeared to be the sum of two or three blended partial 
detection’s. Hence, the MIPS flux was replaced with an up¬ 
per limit, consistently with the “quenched nature” suggested 
by the optical/near-IR colors of the galaxy. 

• Objects #3230 and #3231. These objects are 
blended with a third lower-mass (~ 2 x 10 10 Mq) neigh- 



Figure Al. Comparison between MIPS/24 pm and HST/WFC3 
11-band (F160W) images for objects #427 and #428. Both 
cutouts are 12 arcs in size. The red crosses correspond to 
WFC3/HST positions, while the large open circle shows the MIPS 
FWHM (5"9), centered on the peak of the 24 pm/MIPS emission. 
The open diamonds mark the positions of IRAC sources in the 
field. 


bor (#3232, z p hot ~ 2), not included in the MIPS GOODS 
catalog, which instead seems to be responsible of the major¬ 
ity of the MIPS flux (#24 = 109.6±29.3 pJy). As detailed in 
Appendix m while the object in the middle could be just a 
projected neighbor, #3230 and #3231 are spectroscopically 
confirmed at the same redshift. As also proved by their spec¬ 
tra, #3230 and #3231 seems to be quenched galaxies, with 
MIPS fluxes < 46 pJy and < 62 pJy (sSFR< 10“ lo t/r _1 ), 
respectively. 

• Object #5530. This MIPS source, previously consid¬ 
ered as detected above 3a, was decomposed in three different 
sources, based on the WFC3 positions. A marginal 2a de¬ 
tection is found for the SE component, while the MIPS flux 
of #5530 is below the 2a upper limit. 

• Object #5415. This galaxy was considered as MIPS- 
detected in the automated MIPS-source extraction proce¬ 
dure, based on IRAC priors. However, it is blended with a 
close object (i.e., #5414, at z p hot = 1.96, and not included 
in our sample because of its mass two times below the mass 
limit). The PSF fitting method reveals a secure MIPS detec¬ 
tion (SNR~ 4) for the neighbor, and a marginal detection 
(SNR~ 2.5) for #5415. In this case, we classified #5415 as 
MIPS-d, due to its blue colors (typical of sBzK and blue 
star-forming UVJ), and to the SED fitting results, show¬ 
ing reddening Ay ~ 1, and SFR(SED) comparable with 
SFR(IR±UV), derived from MIPS (cf. Table 0. 

• Object #7923. This galaxy is part of a complex 
blended system, and does not have MIPS counterparts 
within 2 arcsec in the MIPS catalog based on IRAC positions 
(the closest one is 2.4 arcsec away, SE, to this object). The 
MIPS flux decomposition based on WFC3 positions con¬ 
firmed that this galaxies is not detected in MIPS, with a 
SNR< 1, thus in Figure [2] we replaced the 3a upper limit 
with a 2a upper limit . 


APPENDIX B: CLOSE PAIRS, GROUPS AND 
MERGERS 

Since many projected pairs are included in our sample, for 
each galaxy we checked the redshift of all the closest neigh¬ 
bors (dist^ 5 arcsec), so to check if they are also gravita¬ 
tionally bound. As detailed below, some of the bound galaxy 
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Figure A2. Same as figure IaTI for objects #686 and #687. 




Figure A5. Same as figure [ATI for object #1187. 



Figure A3. Same as figure lAll for object #720. 




Figure A6. Same as figure lAll for objects #3230 and #3231. 


pairs in our sample are possibly included in the overdensity 
at z ~ 1.61. In the following we summarize the results for 
the most probable cases, including (i) pairs with at least one 
member with spectroscopic redshift and a companion with 
consistent photometric redshift, (ii) pairs with two photo¬ 
metric redshifts in agreement within A z ~ 0.1, or -just in 
one case- (iii) pairs with visible tidal tails, suggesting inter¬ 
action, although the difference in photometric redshifts is 
lS.z > 0.1. 

• #3231, #3230 are two MIPS-u galaxies with Early- 
Type morphology, at a projected distance of 3.62 arcsec, 
and spectroscopicall y identified at z 3vec = 1.610 and z 3pec = 
1.615, respectively ilMignoli et al.ll2005l . i.e., K20 Survey). 
There is a third bluer galaxy in the middle of this pair, 
with a photometric redshift of z p hot= 2.09, which is the 
closest counterpart for the 24 /rm/MIPS source. Recently, 
iTanaka et al.l d2013| j identified these galaxies, together with 
the X-ray source #3258 in our sample (at z = 1.605, and 
23 arcsec far from #3231), as likely members of an X-ray 
detected group of 8 galaxies, located ~ 5’ (2.5 Mpc) away 
from center of the highest-density region in the overdensity 


at 2 ~ 1.61. Object #3231 has been recognized by the au¬ 
thors as the brightest group galaxy (BGG). It is interesting 
to note that its size is comparable to that of local ETGs 
with similar mass. 

• #686 and #687: The PACS-d #686 is classified as 
BLAGN/QSO-1 in the li terature, and spe ctroscopically con¬ 
firmed at 2 = 1.617 dSzokolv et al.l [20041 . see also Table 1, 
and Section m - The quenched companion #687 has an 
early-type morphology, lies at a projected distance of d~ 
l” 07 and has a “likely” sp ectroscopic redshift of z spec — 
1.609 dVanzella et al.ll2008i f. Based on the criterion used in 
the previous section, these galaxies would be part of the 
highest density peak in the overdensity at 2 = 1.61. How¬ 
ever, if both redshifts are correct, these objects would have 
a comoving radial distance of ~ 12 Mpc from each other, 
hen ce wou ld not be int eracting. 

• iKurk et al.l d2009lf identified a close triplet of galax¬ 
ies, including #986 (z=1.609), indicated as the brightest 
confirmed galaxy member, #1084 (z=1.614), and #880 
(z=1.612). All three objects are quenched/quenching ETGs, 
with the largest distance between them of 20 / /0 (169 kpc) 



Figure A4. Same as figure I Al I for object #880. 


Figure AT. Same as figure lAll for object #5530. 
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Figure A8. Same as figure lAll for object #5415. 



Figure A9. Same as figure lAll for object #7923. 


and with the largest difference in redshift corresponding to 
~ 575 km s- 1 . They lie about 1.5 arcmin (i.e., 760 kpc) 
from the center of the h ighest densi t y reg ion described by 
Castellano et al. ('20071. feurk et all 1200^ 1 concluded that 
the relatively small crossing time of a galaxy in this triplet 
(~ 3 x 10 8 yr), indicates that they are possibly undergoing a 
merging. 

• The MIPS-u #1187 galaxy ( z p hot = 1.66) has two very 
close companions. While the closest object at a projected 
distance of 1.27 arcsecond has z spec = 0.86 (IVanzella et al.l 
120051 . 12 OO 6 I . l2008h . the second one (at a distance of 2.17 arc- 
second) has Zphot = 1-54 (ID:GOODS-MUSIC-15228). This 
source seems to be a quenched galaxy (being also undetected 
at MIPS/24 /im), with red colors, and spheroidal morphol¬ 
ogy (n = 2.33, Retire ~ 1.2 kpc). Moreover, also a bluer, 
fainter galaxy (#1188 in the D07 catalog) at z p hot = 1-57 
lies at a distance of 4.27 arcsecond from #1187. Hence, it 
may be that these three galaxies are gravitationally bound, 
and may even be part of the overdensity at z = 1.61. 

• The highly obscured PACS-d galaxy #5534 with 
Zspec = 1.616 has an asymmetric morphology indicative of 
a merging system. Moreover, it has a bluer neighbor at a 
distance of 3.02 arcsecond, with z p hot =1.59 (ID:GOODS- 
MUSIC-6738). 

• Object #720 ( z p hot = 1.89) shows a tidal feature to¬ 
ward the closest neighbor located at a projected distance of 
~ 2 arcsecond, which becomes even more evident when the 
derived 2D surface brightness best-fit model is subtracted to 
the real galaxy image (e.g., GALFIT residual image). The 
putative companion is a massive galaxy (M* ~ 5x 10 11 Mq) 
with visible spiral ar ms, with an insecure spectroscopic red- 
shift Zspec — 2.49 dPopesso et al.|[2009tl . in disagreement 
with the photometric one, designed as the best solution, at 
z p hot = 1-43 (ID:GOODS-MUSIC-16094). Difficult to say 
whether the two galaxies are a physical pair or not. 


• Objects #427 and #428 represent the most peculiar 
pair in our sample, since it seems to be a double lens system, 
as detailed in Section 15.31 Although #427 is spectroscopi¬ 
cally confirmed at 2 = 1.427 while #428 has a photometric 
redshift z p hot = 1.59, they may represent an interactive pair, 
as suggested by the bridge between them as discussed in Sec¬ 
tion [573] 






























